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ON THE OPTIMUM TRAJECTORY OF A ROCKET * 


BY 
ROBERT R. NEWTON ' 


ABSTRACT 


In this paper, the problem of optimizing some functional of a rocket trajectory is 
examined. It is shown that an optimum trajectory may be found if no more than 
four (in a two-dimensional theory) auxiliary conditions are specified. The trajec- 
tories to give maximum range and to put a satellite into an orbit most efficiently are 
calculated, in a vacuum and including effects of the atmosphere. The effects of air 
for the satellite problem are discussed only qualitatively, since the method of allowing 
for them is quite similar to that used for the maximum range problem. 


1. INTRODUCTION 


Let the thrust of a rocket be some arbitrary function of time, T(é); 
and let ¢;(#), ¢2(t), be two coordinates which give the orientation of the 
line of action of the thrust. Then, in an external force field of known 
properties, the three functions (7, ¢:, ¢:) determine the trajectory 
r(t). We want to choose (T, ¢:, ¢:) so as to make r(t) an optimum, 
in some sense. 

The collection of optimum trajectory problems can be categorized 
in various ways: 


(a) By the function of the rocket. The maximum range problem 
and the orbit problem (that is, putting the rocket at the end of burning 
into a specified orbit, which may or may not intersect the Earth’s 
surface) are two examples. Each of these problems may be phrased 
in either of two ways: We may seek to maximize the range, maximum 


* Most of this work was performed under a contract between Tulane University and 
Redstone Arsenal, as a part of Project TU2-9F. Some revision, mostly editorial in nature, 
was carried out at The Applied Physics Laboratory, under Contract NORD-7386, Navy 
Bureau of Ordnance. 

1Formerly, Department of Physics, Tulane University, New Orleans, La.; presently, 
Applied Physics Laboratory, The Johns Hopkins University, Silver Spring, Md. 
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altitude, or other such property, for a given fuel consumption, or we 
may specify such a property and seek to minimize the fuel consumed. 

(6) By the degrees of freedom allowed. If T(¢) is fixed, and only 
the orientation is varied, we shall call this the restricted problem. If 
both 7 and the orientation angles may be varied, we have what we shall 
call the general problem. Between these two problems lies what we 
shall call the burning-time problem. In this problem, the shape of the 
thrust curve is specified, but not the scale; that is, T(#) = f(at), where 
f is known, but the scale factor a is free. 

(c) By other attendant circumstances. We may assume a flat or 
curved Earth; if curved, it may be rotating or non-rotating. We may 
neglect or include aerodynamic forces. 


Fried and Richardson (1)? have discussed the restricted maximum 
range problem, with both a flat and curved Earth. Fried and Lawden 
(2) have discussed the restricted orbit problem in a vacuum. Several 
writers (3, 4, 5, 6) have discussed the general problem of reaching maxi- 
mum altitude in vertical firing. Bryson and Ross (5) also considered 
non-vertical firing, but with the special assumption that the angle of 
attack of the rocket is always zero. Lawden (6) also considered non- 
vertical firing, but neglected the dependence of aerodynamic forces 
upon the attitude, a dependence which we shall show to be important. 
Hibbs (7) has discussed the best thrust function for horizontal flight, 
while Miele (8) has considered the minimization of aerodynamic heat- 
ing effects. 

We shall discuss the maximum range problem for all categories 
under (0) above, with and without aerodynamic effects. The discus- 
sion will include vertical firing for maximum altitude as a special case; 
our results differ slightly from those of earlier writers for reasons which 
will be discussed later. We shall neglect the Earth’s curvature during 
burning; it will become apparent that the distance travelled during 
burning with optimum conditions is so short that this is quite accurate. 

We shall discuss the orbit problem in detail only ina vacuum. The 
techniques for including aerodynamic effects are so similar to those for 
the maximum range problem that their separate detailed discussion for 
the orbit problem does not seem warranted. Our treatment will differ 
from that of Fried (2) in that we shall allow a period of no thrust to 
occur within the over-all burning period, thus effecting a considerable 
saving. The advantages of a ‘‘coasting’”’ period have been pointed 
out by Lawden (2) and others. 

We tacitly assume throughout this paper that the rocket is in near- 
vacuum conditions for most of its unpowered flight. 

It is possible to “‘over-determine”’ an optimum trajectory problem by 
specifying too many boundary conditions. The nature of the allowable 


2 The boldface numbers in parentheses refer to the references appended to this paper. 
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boundary conditions is not obvious, and we shall discuss this topic in 
the next section. A problem may also fail to have a solution for other 
reasons. We have not found any general existence theorem, but in 
some problems with no solution, or at least with no analytic solution, we 
can demonstrate the non-existence by qualitative arguments. 


2. ALLOWABLE BOUNDARY CONDITIONS 


Over a flat Earth, it is obvious that the optimum trajectory during 
burning should lie in the vertical plane containing the position and the 
velocity vector at burn-out. Accordingly, we use two coordinate axes 
in this plane, £ taken horizontal and 7 vertical. Let ¢ be the angle 
from é to the direction of thrust. Further useful coordinates are 


= tan (9/8), (2.1) 
8. (2.2) 


The forces assumed acting on the rocket are: (a) the jet thrust 7, 
(b) gravity, (c) aerodynamic drag, and (d) aerodynamic lift. If M is 
the instantaneous mass, so that — M is the rate of burning fuel, and if 
p, v, A are air density, rocket speed, and a reference area of the rocket, 


we assume: 


T = — (2.3) 
drag = Cop*A/2, (2.4) 
lift = Czpv*A/2, (2.5) 


with v, constant. Cp and C, are functions of Mach number; Cp is an 
even function, and C, an odd function, of a. 
The equations of motion are: 


Mé = T cos ¢ — 3CpApvé — (2.6) 
Mi = T sin — + Apvé — Mg. 
These are to be solved subject to initial conditions: 
£(0) = 7(0) = 0, &(0) = 4(0) = 0. (2.6a) 


Since we assume that the orientation is controlled, we do not need 


any rotational equations of motion. 
In the general problem, the functions M, 7, and ¢ are not known 


a priori. We can solve Egs. 2.2, 2.3, and 2.6 for these functions; the 
solution for M is of the form: 


M = g, é, ii). (2.7) 
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In this section, we are interested only in the arguments of f, and nct its 
form ; let the arguments be denoted, respectively, by M, x1, %2, «++, Xe. 

The optimum trajectory is that set of functions x,(#), 7 = 1, ---, 6, 
which makes 


M(t) = M(O) + f 


a maximum, for a given M(0), and subject to conditions not yet stated. 
Let ¢, be the burning time on the correct trajectory, and let ¢, be the 
burning time on a varied trajectory, with 6¢, = t; — t. Taking varia- 
tions in Eq. 2.7: 


5M = fubM + fedxi, (2.8) 


using a subscript to denote a formal partial derivative. Since 6M is 
independent of M, Eq. 2.8 can be integrated; since we want 6M (t,) 
to be zero: 


= 6M(t,) = 5M + M(ts)dt, 


x (5 fa Jat + ar 


Now we replace each x, by its original symbol, and proceed, in the usual 
fashion, to eliminate each variation of a derivative, such as 5£, using 
integration by parts. This process gives: 


O = M(tr)dt. + + (fufi — (fi) + 


+exp (f° futoat) exp (- + 


fufi + (f° — + fi — (fu) 


+ similar terms in 7. 


The integrated terms contain no part from ¢ = 0, because of Eqs. 2.6a. 

5£(¢) and 6n(¢) are independent, arbitrary, functions of ¢, except at 
the end points. Hence, the coefficients of 5&(¢) and 6n(¢) must be 
identically zero: 


fe — (fi) + fufit (fd — fi + (fu)? fi — (fu) fe = 9, (2.9) 
fr — fad’ + + (fd — + — (fu) fe = 0.) 


6 
|| 
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Also, the entire integrated part must be zero, but not the individual 
coefficients of 5£(t,), etc., because the end-point variations need not be 
independent. If we let 6& = 65é(t,), etc.: 


+ — (f° + + + 
Chis — (f° + + M (te) dt,= 0. (2.10) 


So far, we have not said whether we are treating a maximum range 
problem or an orbit problem. Therefore, Eqs. 2.9 and 2.10, as well as 
2.6a, hold in all cases. Each of Eqs. 2.9 is fourth order, and the pair 
is a system of eighth order, and, in addition, the first order equation (2.7) 
must be satisfied. The optimum trajectory thus satisfies a differential 
equation system of the ninth order, in three variables M, &, 7, whose 
general solution contains nine arbitrary constants, in addition to the 
unknown time interval ¢,. Thus, altogether ten, and only ten, boundary 
conditions may be specified. 

The initial value of M gives one of these conditions; Eqs. 2.6a give 
four more. Thus there are exactly five boundary conditions yet un- 
specified. We shall find that these five are specified by the conditions 
imposed at the end of burning. 

The range is a function of the conditions at the end of burning, that 
is, at ¢ = ¢,; in a rectangular geometry, for example, in a vacuum: 


R = + (&/g)E[m + Van? + 2gm). (2.11) 


In the maximum range problem, we can fix R, and optimize M. In 
other problems, we fix different functions of 91, Thus, for 
either type of problem, we require s functions of the coordinates and 
velocities at the end of burning to vanish: 


g,(&1, My m) 0, 2; (2.12) 


where the value of s depends upon the problem. 

When we vary the trajectory, &; = & + 6&:, etc., where & is the 
value for the correct trajectory, and = + etc. The change 
in g, is thus 


5g, = (0g,/081) + Edt.) + --- 
+ (8% + = 0, r= 


This set of s equations can be solved for s of the variations 5£,, - - -, di, 
dt,, in terms of the remaining (5 — s) variations,* and the results sub- 
stituted into Eq. 2.10. The (5 — s) variations remaining in Eq. 2.10 
are now arbitrary and independent, so the coefficient of each must van- 
ish. Thus, Eq. 2.10 gives us (5 — s) boundary conditions, and Eqs. 


3 Hence, s cannot exceed 4, since the equations are homogeneous. 
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2.12 give us s more, making 5 in all. This exhausts the allowable 
number of boundary conditions. 

To summarize the general problem, then, we minimize the fuel 
consumed, subject to specifying the following conditions: the initial 
mass, position, and velocity, and no more than four‘ conditions on the 
position and velocity at the end of burning. 

In the restricted problem, we must proceed somewhat differently 
We must first recognize that it does not make good sense to specify 
that the functions g, of Eqs. 2.12 shall vanish. For example, a rocket 
with a given M(t) function might be capable of a range of 6 X 10® ft. 
We could certainly find an orientation function ¢(¢) such that the func- 
tion g = R — 3 X 10° vanishes, that is, such that the range is half of 
the potential range, but this does not tell us anything about how to use 
the given rocket most efficiently. Hence, we shall always phrase the 
restricted problem in terms of maximizing some function of the condi- 
tions at burn-out. 

We can also phrase the general problem in similar terms, by fixing 
the mass M at the end of burning, maximizing one of the g’s, say g1, 
while still fixing ge, ---, g.. The discussion of the general problem is 
altered in only one trivial respect. Equations 2.9 still hold, we get 
(5 — s) conditions from Eq. 2.10, (s — 1) from fixing ge, ---, g,, and 
one from fixing M;. 

Since the restricted problem is a special case of the general problem, 
it cannot allow any more freedom in choosing boundary conditions. 
Also, by inspection of the allowed conditions for the general problem, 
we see that the restricted problem allows all of them. Hence, the 
same boundary conditions are allowed for both problems. 

This conclusion obviously holds for the burning-time problem as well. 

We have assumed a rectangular coordinate system; however, the 
arguments used make use only of the number of coordinates, and the 
orders of their derivatives which occur, and hence do not depend upon 
the coordinate system. If we generalize to three dimensions, we add 
another fourth-order differential equation to Eqs. 2.9 ; we add two terms, 
and hence get two more boundary conditions, in Eq. 2.10; and have one 
more coordinate and velocity to specify initially. 

However, we can now allow s, the number of conditions at burn-out 
in Eqs. 2.12, to be increased to six. This is the only change needed in 
the earlier statement of conditions which may be specified, for any 
coordinate system in three dimensions. 


3. MAXIMUM RANGE PROBLEMS 


The problems of maximum range, whether general, restricted, or 
burning time, are capable of an almost endless amount of variation, and 
we shall not pretend to be exhaustive. Also, except in a few simple 


4 For three-dimensional motion, this is increased to six conditions. See the last paragraph 
of this section. 


Sept., 1958.] OptimuM TRAJECTORY OF A ROCKET 161 


cases, we cannot expect to get an answer in closed form. What we do 
hope to accomplish is to obtain methods of approximation which can 
be applied in a wide variety of cases. 


a. The Restricted Problem in a Vacuum 


This problem has been solved by Fried and Richardson (1) ; we shall 
apply the results, in our notation, to the special case of constant thrust 
and constant gravitation. The central result for this problem is that 
¢(t) is constant, given by 


tan @ = + (3.1) 


This result holds for any form of the function T(#). The values of 
and of course depend upon 7 (¢), as well as upon the value of 
¢, so that Eq. 3.1 is a transcendental equation to be solved for ¢. 

For constant thrust, M is constant, and: 


M = M, + Mi, 


T = — M2,, M< 0. (3.2) 


Letting 
w=1- Bh, (3.3) 


the position and velocity at the end of burning are readily found to be: 


— v, cos ¢ In ps, 


= cos In ws + 1 — (3.4) 
= — v,singdln — g(1 — 
no = [v. sin In ws + 1 — — g(1 — 


If we substitute these into Eq. 3.1, we find® 
(cos 2¢/sin* ¢) = [2(In ws + 1 — ws) /In® ](g/Bv.). (3.5) 


Since cos 2¢ = 1 — 2sin?¢, this is a cubic to be solved for sin ¢. 
Substituting for 8 from Eq. 3.5 into Eqs. 3.4 and then into 2.11, we 
obtain 


gR = (v2 In? (cos ¢/2 sin’ ¢). (3.6) 


Since In ws + 1 — ws in Eq. 3.5 is negative, cos 2¢ must be negative, 
and @ > 45°, as we expect. 
As a numerical example, suppose 
ws = 0.19, t, = 50sec., v, = 7000 ft./sec., 
B = 0.0162 sec.", = — 1.66073. 


5 In? uw, means (In 


(3.7) 


el 
al 
1e 
¥ 
t 
f 
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Substituting into Eq. 3.5 and solving for ¢, we find 


@ = 0.819 622 rad. = 46.96°, 
sin @ = 0.730 888, cos ¢ = 0.682 497, 


using g = 32.2 ft./sec.* exactly. For the range, we get 
R = 366.8 X 10* ft. 


A few other quantities are of interest. v, is the initial acceleration 
due to the jet, about 3.5g, The initial horizontal acceleration is 2.40, 
the vertical is 1.58g. Just before the end of burning, the jet accelera- 
tion is 598 ft./sec.*, or 18.6g. The burnt velocities are § = 7954 
ft./sec., 7, = 6887 ft./sec. The trajectory at burn-out makes an angle 
of 40.9° with the horizontal. This is less than the angle which maxi- 
mizes the range after burning, so in maximizing the range, it is not 
correct to impose as conditions at the end of burning those conditions 
which maximize the range of the corresponding inert projectile. 

It is usually adequate to neglect £, and 7; in Eq. 2.11, yielding: 


gR = 
Using this, 
(cos 2p/sin ¢) = (g/Bv.)(1 — ws)/In ps, (3.9) 


gR = 2v2 In? (cos* ¢/sin ¢), (3.10) 


in place of Eqs. 3.5 and 3.6. Using the numerical values of Eq. 3.7, we 
get ¢ = 47.95° instead of 46.96°. If we used the incorrect ¢, we should 
get an actual range of 366.6 X 10‘ ft., a loss of less than one part in a 
thousand. 


b. The General Problem in a Vacuum 


As is well known, the equations for the general problem in a vacuum, 
that is Eqs. 2.9 and 2.10, have no analytic solution. We can easily see 
this by looking first at the burning-time problem, with the aid of Eqs. 
3.5 and 3.6. In Eq. 3.6, for a constant value of v, In uw, the greatest 
range is attained when ¢ = 45°. This does not happen because there 
is a relative maximum of R at 45°; instead, 45° is the smallest value of 
¢@ with physical significance, and R decreases as ¢ increases. From 
Eq. 3.5, when ¢ is 45°, 8 is infinite, and ¢; is zero. Thus there is no 
analytic solution ; the closer we can come to instantaneous burning, the 
better. 

If t, is zero, &, and are zero, while £, and 7 are not. The range 
is then 


R = 2&s%/g. (3.11) 


U. 
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If we reach the condition of instantaneous burning as the limit of any 
other burning curve than a constant, it does not affect the results, 
because the velocity depends only upon », and the ratio of initial to 
final mass. Therefore, all modes of instantaneous burning give the 
same range. 

If we take a more inclusive approach, we shall find that there is an 
optimum burning curve, even ina vacuum. Suppose we have a rocket 
with constant burning rate, having a long burning time. We can 
improve the range by shortening ¢,, but not as rapidly as Eqs. 3.5 and 
3.6 indicate: As we shorten ¢,, we must increase the burning rate, the 
motor pressure, and the motor weight. Eventually, the loss from 
increased motor weight out-weighs the gain from a shorter ¢,, so there 
is an optimum ¢,. Further, the optimum shape of the burning rate 
curve is probably nearly constant ; if it were not, the motor weight would 
be determined by the maximum burning rate, and most of the time 
we should therefore be carrying excess weight, and losing range as a 
result. 

The loss of range from having a non-zero burning time is readily 
estimated from Eqs. 3.5 and 3.6, provided the parameter $v./g, the 
initial ratio of thrust to weight, is not too small. Assuming that 6v,/g 
is large, we can expand Eggs. 3.5 and 3.6 in powers of ¢ — (2/4), keep 
only terms through the first power, and eliminate ¢, obtaining: 


gR = 02 In? — (g/Bv.)[v2(us — 1 — In ws)/In* (3.12) 


This is quite accurate even for Bv./g as small as 3. We see that the loss 
in range varies inversely as 6v./g; for our example, this loss is 12.6 


per cent. 


c. The Restricted Problem in Air 

™ The first question which one would probably ask about the problem 
in air is the accuracy to which the vacuum problem approximates it. 
One can estimate this accuracy by estimating the aerodynamic forces, 
using velocities calculated from the vacuum theory ; the most important 
aerodynamic force for present purposes is the drag at zero angle of 


attack. 
As the rocket climbs through the air, its velocity increases, while the 


air density decreases exponentially. We shall write the air density p as 
p = psxf(n), (3.13) 


where psx is the density at sea level. For convenience, we shall 
assume : 


f(n) = exp(— an), (3.132) 


. 
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where a = 3.34 X 10-5 ft.-'; this is accurate up to 35,000 ft. Since the 
air density decreases so rapidly, the drag force should reach a maximum, 
which we wish to estimate roughly. 

Assuming a straight trajectory at an angle ® with the horizontal, 
traversed at constant acceleration G, and a drag proportional to some 
power, a say, of velocity, we want the maximum of f* exp(— aG# sin @/2). 
This maximum is found to occur always at the same altitude, a/2a, 
independent of G. The velocity at this altitude, which is about 20,000 
ft., is almost surely well into the supersonic region ; in this region, a = 1 
with reasonable accuracy over a considerable velocity range. Writing 
the drag as gpvf(n), f(n) is about 0.6 when drag is a maximum, and the 
maximum is about 0.6gp(G/a sin @)'*. We now need some estimates 
for gp. 

For a rocket of cross-sectional area A, and a drag coefficient of 0.2 
at Mach number 4, we get gp = 30A in f.p.s. units. This must be 
compared with the mass in order to be significant ; for an over-all density 
of 100 Ib./ft.* for the rocket, and a length-to-weight ratio of about 10, 
the mass M is about 10*A*/*, Thus: 


gp 0.3 
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For the ratio, maximum drag divided by weight, we correspondingly 
obtain : 


using f.p.s. units for G and M. For rockets capable of reaching high 
altitudes, this ratio is not likely to exceed unity, hence the ratio of 
drag to thrust for such rockets is expected to be much less than unity. 
From the preceding discussion, we can reach two conclusions: First, 
the vacuum theory gives a reasonable approximation for a wide variety 
of cases. Second, aerodynamic forces are probably most important at 
supersonic velocities. In accordance with the latter conclusion, we 
shall find it convenient to write: 


gof(n)v(1 + 
gxf(n)va, 


drag = 


lift 


with f() defined by (3.13), and with gp and gz being slowly varying 
functions of velocity. The constant kp is a measure of induced drag; 


kp is approximately equal to 1 + (g:/gp). For numerical examples, 
we shall use 


= 0.0162 Mo Ib./sec., 


gu =S5gn, kno = 6. (3.15) 


These values were chosen for convenience, and seem to us to make 
aerodynamic forces about as large as they are likely to be for the class 


A 
: 
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of problems we are considering. In any actual case in which aero- 
dynamic data are available, one should either use actual values, or use 
the best fit to the form of Eqs. 3.14 at the velocities where drag is 
expected to be largest. 

We shall continue, in this section, to neglect drag after the end of 
burning, and to use Eq. 2.11 for the range. Unless the drag after 
burning is quite high, the conditions which maximize the range in 
vacuum should be close to those which maximize the actual range. 

If the vacuum theory furnishes a reasonable approximation, what 
we must do is to calculate differences from the vacuum theory. This 
suggests that we transform to the coordinate system of Fig. 1, with an 


+ 


Fic. 1. Coordinates used for the maximum range problem in air. 


x-axis rotated from the é-axis by an angle ®. For #, we use the con- 
stant value of ¢ furnished by the vacuum theory. We have chosen the 
y-axis in the direction shown so as to make y positive for the vacuum 
solution. Let x be the angle between the rocket axis, which is the 
direction of the jet thrust 7, and the x-axis, and let a be the angle of 
attack. All angles are positive if they have the senses of rotation 
shown in the figure. If 6 is the angle between x and the velocity 2, 
so that 


6 = = y/v, (3.16) 
the angles 6, x, and a are related by 
6=x+a. (3.17) 


a 

ae 
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If the (x, y) coordinate system, the equations of motion (2.6) be- 
come, upon using Eqs. 3.14 for drag and lift, and after dividing by the 
mass M, and letting G = T/M 


= Gcosx — gsin® — M-f(n)gn(1 + kpa*)t + Mf (n) gray, (3.18) 
j = Gsin x + gcos® — M-f(n)gn(1 + kva*)y — Mf (n) 


the transformation from x, y to é, n, being: 


+ ysin ®, (3.19) 
n = xsin ® — ycos ®. 


We expect y and all its derivatives to be small, and therefore also 
expect the angles 6, x, and a to be small. Therefore, we take as a zero- 
order approximation to the trajectory 


= G — gsin® — (3.20) 


in which 7 is x» sin ®@. We have written gp(zo) to emphasize that gp is 
a function of velocity, and that the velocity is approximated by Zo. 
Equation 3.20 cannot be solved by quadratures, but it can readily be 
integrated numerically, or solved by a perturbation procedure. In the 
rest of this portion of the paper, we assume that this has been done, so 
that x(t) is a known function of time. 

For brevity, let us use the notations: 


D = go(%o)f(0)/M, L = (3.21) 
N=L+4+D, D’ = dD/dno; 
all of these quantities are known functions of time through our knowl- 
edge of xo. Let us further define three quantities yo, «1, and y; by means 
of the differential equations 
Got Nyo=Gx+g cos $+Liox, 
#:+Dé,= — (Gx?/2) +40 cos ®D’ ao + LO — — Dk viva’, (3.22) 
— Dk pyre’, 
with 6) = #o/%o. These equations can be solved successively by quad- 


ratures. By direct substitution, remembering Eq. 3.17, we can verify 
that 


x = xo(t) + x(t), = yo(t) + (3.23) 


order of small quantities. We assume that gp(v) can always be re- 


are the solutions of the equations of motion (3.18), through the first 
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placed by gp(#o); this should be highly accurate, since gp varies slowly 
with v, and the y component of velocity hardly alters the speed. 

By repeated integrations in the first of Eqs. 3.22, we can obtain 
successively #o(¢) and yo(t). Then, by substitution into the other equa- 
tions, further integrations yield x,(¢) and y:(#). By systematic reversal 


of the order of integration in all multiple integrations, we can write 
all results in terms of single integrals, if we make use of the functions 


Iy(t) = exp [- Iyi(t) = In(t:)dt,, (3.24) 
and similarly defined functions Jp(#) and Jp'(t). The results are: 


yo = f [Ini (t) — (ts) + gcos® + Léox (3.26) 


X [D’ cos Syo%o + LOoyo — Lox — — (3.27) 


- 


x,= 

X [D’ cos + — — 3Gx* — (3.28) 
= — f° (3.29) 
yt(t) — In*(ts) (ts) "Dk (3.30) 


Dummy variables of integration have been omitted where there seems 


little danger of ambiguity. 
By taking variations in Eq. 2.11 


where 
gt; = no + (me? + (3.31) 


t, being the time of flight after burning. It is clear that 


= 6x1, cos + (dyon + dy) sin ®, (3.32) 
dns = 6x1, sin — + cos 


- 
> 


168 Rosert R. NEwTon 


since x» does not involve x, and hence 6x» is zero. Then: 


= 0 = {cos® + + sin + 
+ {sin® — &5( + cos {yon + + ty + 


We adopt further abbreviations: 


A = (no? + 2g.) cos ® + &, sin , (3.33) 
S = (a? + 2gm)'? sin @ — £, cos ®, 


in terms of which 
A (5x15 + + + dyin + + = 0. 


A and S are not known, because the terminal position and velocity 


depend upon x. 

We get the variations dyo,, etc., from Eqs. 3.25 to 3.30 by a succes- 
sion of simple steps: Replace ¢ by ¢s, allowing us to replace ¢; by ¢ without 
ambiguity, then take the variations of the integrals. Since the limits 
are fixed, we can take variations under the integral sign ; quantities in 
the integrands which are to be varied are x, Yo, Yo, 90, anda. The only 
independent variation is 6x; for the others, we take 


= Iv(t) + Lav) dxdt, 


= [Iwi(t) — (ts) + Lao)ixdt, 


500 = Sa = 500 — dx, 


it being understood that the last two are to be expressed in terms of 5x, 


using the integral for dy. 
When we carry out all of these steps, we finally arrive at a single 


integral which is equal to zero, and whose integrand is the product of 5x 
with a complicated function involving xo, yo, x, etc. Since dx is arbi- 
trary, the function which multiplies 5x must equal zero. Using still 


another abbreviation : 
In (to, t) = — Int (t) + tI w(t) (3.34) 
and a similar abbreviation 9 p(t,, t), this function is 
— Ad p(to, t)(Gx + Lio) + + Lao) + t) 
X — x) + (te, — x) + Iw) + Lao) 
x f° (20 — x) — 2Dko(O» x)] 
+ [Iw'(t:) — p(t, t1)D’ cos = 0. 


UJ. F. I. 
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S is small compared with A, and y is small compared with z. Hence, 
the fourth term can be neglected in comparison with the third. The 
function x(#) must now be chosen so as to satisfy this equation. 

This equation is a complicated integral equation for x. Not only 
does x occur explicitly, but also implicitly, since the functions Yo, 0, 9, 
and gp, as well as the constants A and S, depend upon the choice of x. 
The equation is best solved by successive approximations ; a convenient 
form of the equation for this purpose is® 


In9n(G + + Lao)-x 
= + f (ti) — In'(t)]9 t1)D’ cos 
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+ — L)(G + + f 
X[L(x + 2a) — 2DRpa (3.35) 


None of the terms on the right is really independent of x, however, the 
first two depend upon x but slightly. 

Hence, we have the following steps in calculating an approximate 
solution x(t) for the best orientation of the axis: 


1. Find the best value of ¢ in a vacuum from Eq. 3.5, or an equiva- 
lent equation if J is not constant. Call this value ®. 

2. Solve Eq. 3.20. Using the solution x, tabulate the functions 
D, L, D’, In, Ix’, Ip, Ip', defined in Eqs. 3.21 and 3.24. 

3. Assuming that y is identically zero, and that x equals xo, calculate 
A, S, gn(ts, t), and gp(t, t) from Eqs. 3.33 and 3.34. 

4. Neglecting the last two terms in the right member of Eq. 3.35, 
solve that equation to get a first approximation to x, say xo. 

5. Use xo for x in Eqs. 3.25 and 3.26 to calculate yo and yo, and then 
calculate a from a = (%o/%0) — xo. 

6. Using the values from Step 5, calculate the last two terms on the 
right of Eq. 3.35, and combine with the estimate already made of the 
first two terms to get a new approximation for x, say x1. 


We have carried this process out, with the numerical values used in 
our previous numerical examples. The results for xo(¢) and x:(¢) are 
shown in Fig. 2. We see that x, differs considerably from xo for small ¢. 
Such behavior would normally lead us to think that the method of suc- 
cessive approximations used in calculating x; does not converge rapidly, 
and hence that x: itself is a poor approximation. However, the early 
values of x are determined to a considerable extent by the values of x 
late in burning, when xo and x, agree rather well, and x: seems to be 
actually a good approximation to the solution of Eq. 3.35. 

® As we shall see in a minute, x can be replaced by tan x in the left member of Eq. 3.35, 
but not in the right member. 


170 Rosert R. NEwTon {J. F. I. 


However, there is still a possible difficulty. The early values of x: 
may be large enough that the use of x for sin x may not be justified. 
This difficulty disappears if we note that, when ¢ and hence 4» are small, 
the left member of Eq. 3.35 is approximately 


(Gx/G). 


If we trace through the derivation, we note that Gx, in the numerator, 
came from differentiating G cos x, yielding — Gsin x, then approxi- 
mating this by — Gy, and finally changing the sign. The G in the 
0.6 
0.5 
0.4 


0.3 


Fic. 2. An optimum orientation function for the maximum range problem. xo is the zero- 
order approximation, x: is the tangent of the angle given by the first-order approximation. 


denominator came from differentiating G sin x, yielding G cos x, then 
replacing cos x by unity. Hence, x in the left member of Eq. 3.35, 
but not elsewhere, can be replaced by tan x. Thus, as a final step 7 
to follow the steps previously enumerated : 


7. For the optimum x(t), use tan! x:. 


Equation 3.35 equates the change in a loss, the left member, to the 
sum of the changes in several gains, as we try to improve the range by 
pointing the axis up more steeply than the constant angle @. The 
gains are: 


xX, 
0.2 
0.1 
2 t, sec. | 
lO 20 . 30 
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1. Lift, which is upward and tends to increase the range. 
2. Increased altitude, which decreases the drag. 
3. Smaller a, which decreases the induced drag. 


For our example, the effects of increased altitude and of induced drag 
are about equal (xo primarily depends upon the altitude effect, and 
x0 — x1 upon the induced drag effect). If the lift is large, so is the 
induced drag, and induced drag gives the dominant effect of these three. 
The loss is a very severe one; as we point the axis away from the desired 
direction of travel, we acquire momentum in an undesired direction. 
This momentum must be killed by adding yet more momentum in the 
equally undesired opposite direction. Hence, we cannot tolerate point- 
ing the rocket off the desired heading very long. 

In summary, the optimum trajectory in air is determined mostly 
by allowing a small loss in the form of undesired momentum to be 
balanced by the advantages of a small angle of attack, and a steep climb 
of short duration, in order to decrease air density. The optimum 
trajectory, then, tends to be characterized by small angles of attack. 

We took the function x(¢), calculated as above, and used it in a 
numerical integration of the equations of motion. We also integrated 
the equations of motion using the constant tilt angle furnished by the 
vacuum solution, that is, with x(#) = 0. Using x = 0, the resultant 
range was 76.25 per cent of the range in vacuum; the other x(¢) func- 
tion gave 77.56 per cent of the range in vacuum, a gain of 1.31 per cent 
of the vacuum range, or 1.69 per cent of the actual range. 

Such a small gain is to be expected. By using the optimum x(é), 
we do not affect the gravity loss, only the drag loss, which is 23.75 
per cent of the vacuum range if we use the vacuum theory. Even if we 
could fire vertically, we would only expect to reduce the drag loss by a 
factor of about v2 (by decreasing the distance travelled in air), saving 
about 7 per cent of the vacuum range. Since we cannot hold the axis 
vertical for very long, the best gain we can expect is a small fraction of 
7 per cent. 

Because the gain which can be expected over the results given by 
vacuum theory is typically so small, we feel that it is not important to 
have an accurate theory in allowing for aerodynamic effects. This 
constitutes the ultimate justification for some of the assumptions we 
have made, which would be considered quite crude for some purposes. 

Figure 3 shows the shapes of the trajectories using the optimum x (#) 
and using x = 0, the vacuum theory result. The optimum trajectory 
has only a slight curvature, a result which we expect to hold unless the 
initial thrust-to-weight ratio is very nearly unity. The x-axis is also 
shown in Fig. 3 for comparison. 

In this portion of Sec. 3, we have considered the thrust to be a known 
function of time, and varied only the orientation. In the remaining 
parts of this section, we consider both the thrust and the orientation to 


A 
i 
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be variable. However, the results of this section allow us to separate, 
to satisfactory accuracy, the variations in thrust and orientation : Since 
the optimum trajectory is nearly straight, we can study the thrust 
variations by assuming the rocket to be constrained to a straight line, 
and optimize the range under this constraint. Using the thrust func- 
tion so determined, we find the best orientation function for this thrust 
program. If we wish, we can repeat the process, using the newly 
determined curve as the curve of constraint, and so on, but it seems 


12 x 


10 


4 6 8 


10* 


Fic. 3. Shape of an optimum trajectory for the maximum range problem. Curve A is 
the optimum trajectory; Curve B is the trajectory which would be followed in air by a rocket 
using the orientation function which is optimum in a vacuum. The cross-marks on the curves 
are at time intervals of 10 sec. The x-axis is shown for comparison. 


unlikely that this repetition will be warranted in most cases, particu- 
larly since we shall find that the optimum initial thrust is many times 
the initial weight. 

Through the first step of this process, the problem of studying 
the thrust function to give maximum range is identical with the 
problem of attaining maximum summit altitude in vertical firing, 
provided that we use the appropriate component of gravity and the 
appropriate dependence of air density upon position. If we use (3.13a) 
for the variation of density, the latter means simply changing the con- 
stant a in the exponent. Hence, in the remaining maximum range 
problems, we shall speak as though we were studying vertical firing. 


| 
8 | 
| 
B 
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d. The Burning-Time Problem in Air 


In this problem, we start with a given shape of the thrust-time curve, 
and vary the burning time while keeping the total fuel consumption 
constant. 

If the total distance travelled by the rocket is large compared with 
the burning distance and with the effective thickness of the atmosphere, 
the losses can be considered as applied impulsively, and maximization 
of distance is equivalent to minimizing the non-conservative losses in 
velocity. These losses are two: (1) gravity loss during burning, which 
is simply gt,, and (2) the total velocity loss to drag. 

To study the latter loss, let ® denote the drag force, and let: 


M/My,, 
D = Cv* exp (— as), (FoF) 


where s is distance travelled. The velocity loss is the time integral of 
2/M over the time interval for which drag is significant ; because the 
density decreases so rapidly, we can make the upper limit infinity with- 
out appreciable error. 

We therefore need to choose the burning time so as to minimize 


gt + (C/M») f * exp (— (3.37) 


For a given shape of the thrust-time curve, the entire curve can be 
identified by a single parameter 8. f, is then a function of 6; v and yu 
can be considered as formal functions of s and 8. We want the deriva- 
tive of expression (3.37) with respect to 6 to vanish’: 


g(dts/d8) + (C/Mz) exp (— as){C(a 
= 9. (3.38) 


(dv/08) and (du/08) are to be evaluated at constant s; after the end 
of burning, 04/08 = 0. It will usually be accurate enough to use values 
of v and uw from vacuum theory ; if this is done, Eq. 3.38 is readily solved 
numerically. 

In order to proceed further, we need to assume an explicit thrust 
function. We shall again assume constant thrust, and shall let the 
parameter 8 be the 6 of Sec. 3.4. We shall also use a = 1, re- 
placing C by the gp of Sec. 3.c. v and s are the same as 7 and 7 
in Eqs. 3.4, if we let sin @ = 1, use the component of g along s, and re- 
place yu, by w. Substituting into Eq. 3.38, and transforming the vari- 
able of integration to p, 


7 As we shall show in Sect. 3.e, C is discontinuous at the end of burning. This introduces 
another term into the left member of Eq. 3.38, which we shall ignore here. 
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1 

1 — = (gov-/Mog) exp (— as) 

X + 1 — w — (g/Bv.)(1 — (3.39) 
We can write as in the form 

as = (av.?/g)(g/Bv.)[u ln + 1 — — (g/2p0.)(1 — 

If the trajectory is a straight line in an arbitrary direction, a and g 
depend individually upon the direction, but their ratio does not. For 
given us, then, we can evaluate the integral in (3.39) as a function of 


the parameter 6v./g, and use the results either for vertical firing or in 


0.5 
0.4 
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Fic. 4. The integral used in calculating the optimum constant thrust, as a function 
of the initial thrust-to-weight ratio (6v,/g). 


firing for maximum range. Figure 4 shows how the integral varies with 
this parameter for values from unity to 16. Values less than unity are 
not physically significant ; for values greater than about four, for which 
our theory is principally intended to apply, the integral is almost a 
linear function of (8v,/g), the initial thrust-to-weight ratio. We have 
continued to use uw, = 0.19 and v, = 7000 ft./sec. 

For firing at 45° with the horizontal, which gives very nearly 
maximum range, and gp = 0.0162 Mo, as before, the integral must 
equal 0.163. From Fig. 4, 6v./g = 8.35, requiring that 6, equal 
(8.35) (32.2/v2) = 190 ft./sec.*, whence 8 = (190/7000) = 0.0271 sec.-, 
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and ¢, = 29.9 sec. For vertical firing, the integral must equal 0.230, 
making 6v, = 333 ft./sec.*, 8 = 0.0475 sec.—', and t = 17.0 sec. 

It seems plausible that these results give the distance travelled dur- 
ing burning more accurately than the time itself, because the variation 
of the integral with 6v./g is dominated by the appearance of that param- 
eter in exp (— as). In the last example, the altitude at burn-out is 
68,000 ft. If we wish to correct the time to reach this altitude because 
of velocity lost to drag, making a very rough calculation of the correc- 
tion, we change f, to 18.4 sec. and 8 = 0.0441 sec. 

As a check, we integrated the equations of motion numerically for 
8B = 0.0441 sec—. To facilitate comparison with the results of the next 
part, where we calculate the optimum thrust function, we calculated 
the velocity at an altitude of 108,650 ft., and found it to be 9765.1 ft./sec. 
Since the maximum obtainable with variable thrust, but with other 
parameters unchanged, is shown to be 9970.6 ft./sec. at this altitude, 
we must be quite close to the maximum obtainable with constant thrust. 

The method of this section can readily be extended to include any 
variation in structural weight with acceleration or burning time; all 
that is needed is to differentiate M» in expression (3.37), as well as the 
other quantities which depend upon the acceleration. Investigating 
the dependence of M> upon 8 is beyond the scope of this paper. 


e. The General Problem in Air 


As we pointed out earlier, the problem of determining the thrust 
function to give maximum range is almost identical with the optimum 
thrust problem in vertical firing, provided that we use appropriate 
values for g and for the constant in the air density variation. The 
vertical firing problem has been discussed by earlier writers (3, 4, 5, 6); 
we shall first summarize briefly their results. 

The problem is now one-dimensional, so we now have to deal with 
only one of Eqs. 2.9, we have two conditions at the beginning of burning, 
namely 7 = 9 = 0, and will have three boundary conditions imposed 
at the end of burning. However, the differential equation degenerates 
to second-order, so that there are two boundary conditions too many. 
Thus, there is no analytic solution to the problem. 

That the differential equation is only second-order can be discovered 
by substituting into the general form of Eqs. 2.9, or perhaps more simply 
by deriving the Euler-Lagrange equation for the problem directly. 
Writing D(n, 4) for the drag force, the equation of motion is simply 


This integrates immediately to: 


) 

— Mo, = + g) + 

i 
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Since the integral involves only n and 4, the Euler-Lagrange equation is 


second-order ; it is in fact 
v.(9D/dn) — — D— (4 + + (D/v.)] = 0. (3.41) 


This depends only upon the form of the drag function, and not upon the 
size of the drag coefficient. 

Since there are too many boundary conditions, there must be dis- 
continuities in 9, unless the surplus boundary conditions are auto- 
matically satisfied by the solution of the problem formed by omitting 
them. It has been shown that 4 changes discontinuously from 0 to 
some positive value, jo say, at? = 0. This naturally suggests the use 
of a booster, and Leitmann (4) has investigated the modified problem 
posed by having the booster drop off. It has been shown that there is 
no discontinuity in 4 at the end of burning, but the proof depends upon 
assuming no change in the drag coefficient at the ead of burning. 
However, there are at least three reasons why the drag coefficient 
should decrease after burning: 


1. The payload may be detached from the motor. This may pro- 
duce a large change. 

2. Cessation of the jet may alter the air flow around the rocket. 
This is probably a small effect. 

3. Cessation of the jet decreases the drag coefficient by the ratio of 
the base pressure to the free stream dynamic pressure. 


The last effect can be appreciable. Using M for Mach number 
momentarily, the ratio of static to dynamic pressure is (47MM/*)-; 
for a diatomic gas, this is (0.7 M*)-. The fractional change in Cp is 
1/(0.7 M*Cp) = 1/0.56 M, using our earlier assumption about the 
behavior of Cp. At Mach 5, the change is one part in 2.8, at Mach 10, 
one part in 5.6. 

We must now formulate the problem in a way which allows a dis- 
continuity in 9 at the end of burning. Equation 3.40 still correctly 
describes the mass variation, and correctly charges to total fuel con- 
sumption the fuel required to produce the initial discontinuity 7 in 9. 
Equation 3.41 will still be obeyed by 7 until time ¢,; at this time, let 
7 = 2, 7 = m2, and M = M,. Then we shall allow 4 to change dis- 
continuously to the value 4, and M to the value M,, without changing 
n. There is the obvious requirement that 4, > #2. Let the drag func- 
tion during burning be denoted by 9~(n, 4), and after burning by 
Dt(n, 9). 

We shall speak in terms of fixing the summit altitude, and mini- 
mizing the fuel consumption. To impose the condition of fixing the 
summit altitude, we follow Leitmann (4), and imagine the equation 
of motion after burning to be integrated, and the result expressed as 


4 
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9 = g(n). When 4 is zero, 7 has a fixed value which is the summit 
altitude. Evaluating at the beginning of unpowered flight: 


= g(n2) (3.42) 


is the condition to be imposed. 
When we take the variation of M;, the quantities to be varied are 


ts, and as well as n(¢) and for allt Since we already 


_ know the differential equation of the trajectory, we need consider only 


variationsin the end conditions. Proceeding in the usual fashion, we get 


+ [gM2 + D-(n2, 
+ + Jing = 0, (3.43) 


from 6M, = 0. From taking variations in Eq. 3.42, 
— g'(n2)[dn2 + = 0. (3.44) 


Upon eliminating 5t, between Eqs. 3.43 and 3.44, 5%, and dn, become 
independent, so their coefficients vanish, yielding 


+ gM. + D-(n2, = 9, (3.45) 
fol + ] — — D-(n2, = 0. (3.46) 


These, along with (0) = 0, are conditions to be imposed upon the 
solution of Eq. 3.41; this is the correct number of conditions, since the 


end time ¢, is unknown. 
The solution of the system is straightforward in principle. Leit- 


mann (4) has shown that the solution of 
iL(@D-/d4) +2. D-] — gM D- = 0 (3.47) 


satisfies Eq. 3.41, as well as (3.46). We therefore have only to integrate 
(3.47), with (0) = 0, until (3.45) is satisfied. Since %, can be calcu- 
lated from (3.42) once 72 is known, we have the solution to the problem. 

It is useful to rewrite the boundary conditions in a way which in- 
volves the discontinuity. The acceleration immediately after burning 
can be written as 7g’(n2), hence the equation of motion immediately 
after burning is 


g'(n2) + gM, + Dt(n2, = 0. 
Eliminating g’(n2) between this and Eq. 3.45, we get 
+ D- (m2, = Lg + Mo (3.48) 


“4 
O- 
n 
iS 
: 


178 Rosert R. NEwTon {J. F. 1. 


If ®- and 9+ are positive and increase monotonically to infinity as 
7 — ©, it is not difficult to show the following, regarding 7, as the 


unknown in Eq. 3.48: 


1. If D+ = D-, Hm» = 7H is the only root, and is double. 

2. If D+ > D- for the same values of the arguments, the only 
possible roots are with 4, < #2, which violates the constraints. Hence, 
there is no terminal discontinuity in 4 for this case, which is unlikely 
to occur. 

3. If D+ < D~ for the same values of the arguments, which is the 
most probable case, there are two roots; one and only one has 7% > #2, 
and is the desired root. 


Continuing with the numerical values we have used before, we 
solved the system of equations numerically for vertical firing, first 
under the assumption that D+ = D-. The initial velocity %. was 
3,008.7 ft./sec. The acceleration was 105.4 ft./sec.* immediately after 
the initial impulse; at burn-out, which occurred at 20.63 sec., 7 = 946.4 
ft./sec.*, 7 = 9970.6 ft./sec. and » = 108,650 ft. This is the case to be 
compared with the optimum constant thrust, which gave 4 = 9765.1 
ft./sec. at the same altitude, a loss of 2 per cent in velocity. Consider- 
ing the higher structural weights which would be required with the 
variable thrust program, the difference would be even less. 

To investigate the effect of changing the drag function, we assumed 
that the drag coefficient changed by 14 per cent, which happened to be 
convenient for numerical reasons. The trajectory was the same until 
the discontinuity occurred, which was at 17.1 sec., with n: = 78,450 ft., 
ne = 7363.4 ft./sec., #» = 10,177.8 ft./sec. At 108,650 ft., 9 was 
9979.5 ft./sec., a gain of about one part in a thousand. We expect 
the change in drag coefficient to be greater than 14 per cent at 7360 
ft./sec., even with no detachment of the payload, but even so, the gain 
from using an impulse at the end of burning will not be large unless 
there is detachment. 

One further interesting property of the optimum solution should be 
pointed out. In Eq. 3.47, if D- is proportional to 4%, 9D-/09 = aD-/9, 
and the equation becomes 


(gM/D-) = a — 1 + (9/2). 


Since 4 is initially rather large, this says that drag and weight are always 
roughly equal for the optimum. This result seems plausible for the 
constant thrust case, also, since gravity loss and drag loss change in 
opposite directions as we change the thrust, but it is difficult to estab- 
lish it rigorously. 

The fact that the optimum variable thrust function gives only a 
slight improvement over the optimum constant thrust shows that, in 


< 
| 
| 
| 


Sept., 1958.] OptimuM TRAJECTORY OF A ROCKET 179 


this problem as well as in the restricted maximum range problem, the 
optimum is very flat ; a wide departure from optimum conditions causes 
a loss of only one or two per cent. This is perhaps our most important 
conclusion, and it has two important corollaries. First, it is not neces- 
sary to use a highly accurate theory in calculating an optimum trajec- 
tory. Second, it is not necessary that a rocket follow the optimum 


trajectory closely. 
4. ORBIT PROBLEMS 


In this section, we turn from the maximum range problem to the 
more difficult orbit problem. There is, of course, a wide variety of 
possible orbits, and we cannot try to treat them all in detail. We 
choose for treatment the typical one of requiring that the rocket attain 
a horizontal velocity of specified magnitude at a specified altitude. 

The orbit problem is more difficult than the maximum range problem 
because it poses a variational problem with auxiliary conditions, requir- 
ing the use of Lagrange multipliers. We shall not try to solve the 
problem in air, but shall make a few comments about the effect of air, 
after discussing vacuum problems. We shall also discuss some of the 


effects of a round Earth. 


a. The General Problem in a Vaccum 


We require that the rocket, at the end of burning, have zero vertical 
velocity, and specified values of altitude, say ,, and of horizontal 
velocity, say §,. Thus, we have three functions of the type used in 


Eq. 2.12, 


g2 = & — & = 0, (4.1) 
gs =m — = 0. 


We wish to minimize the fuel consumption, subject to these conditions. 

It is not difficult to set up the differential Eqs. 2.9 for the optimum 
trajectory, along with the requisite boundary conditions, and it is not 
difficult to show that the system of equations has no analytic solution. 
Instantaneous burning is obviously not the solution, either, because 
with it, 7, equals zero, rather than »,. 

The trouble lies with conditions (4.1), which do not properly reflect 
our desires. We want the rocket to satisfy these conditions at some 
point, but not necessarily at the end of burning. In the case of a flat 
Earth, we can always find a ballistic trajectory, starting from some 
point on the surface, which satisfies the required conditions at its sum- 
mit. This trajectory has initial horizontal velocity ¢,, and initial 


vertical velocity ¥2gn,. If we use instantaneous burning with these 
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values of burnt velocity, we have a more efficient rocket than any which 
burns all the way to the altitude »,. 

Over a round Earth, we cannot use this procedure, because the de- 
sired orbit may not be one which intersects the Earth’s surface. In- 
stead, we use the well-known idea of transfer orbits. The rocket, at 
the Earth's surface, is initially in a circular orbit, because of the Earth’s 
rotation. In order to establish it in a different orbit, we launch it, 
using instantaneous burning, along an intermediate or transfer orbit, | 
which is tangent to both the initial and the desired final orbit. When it 
arrives at the point of tangency with the desired final orbit, we give it 
the necessary increment of velocity, again using instantaneous burning, 
to put it into the final, rather than the intermediate, orbit. 

Even if we neglect the Earth’s rotation, so that the rocket is not, 
properly speaking, in any orbit to start with, it is still not difficult to 
show that the intermediate orbit should be tangent to the Earth’s 
surface and to the final orbit. The losses which result from using an 
orbit which is not tangent to the Earth’s surface have been discussed 


by Ehricke (9). 


b. The Restricted Problem in a Vacuum 


We shall work only in rectangular coordinates, even if we are discuss- 
ing a rocket which is destined to become a satellite. We shall discuss 
the validity of this in part c of this section. 

As we said in Sec. 2, the restricted problem mustal ways be phrased 
in terms of maximizing some function of &, ---, 9s, hence we cannot 
use Eqs. 4.1 as they stand, or as modified by the conclusions of part a. 
However, forgetting the conclusions of part a for the moment, we can 
use the first two of conditions (4.1), and require that 9, be a maximum. 
For any other orbit problem, we can still find one function to maximize, 
while holding two others fixed. 

For arbitrary jet acceleration G(t) and arbitrary orientation ¢(t), 
the coordinates and velocities at the end of burning are: 


= oat, = (t, — t)G cos dt, 


(4.2) 


in = ["Gsinddt gh, m= 0G sin gat — 
0 


When we vary ¢, 4 is to be zero because 7 is to be a maximum; 6é, 
and 4%, are to be zero because £, and 7 are fixed. Using the standard 
procedure in the calculus of variations, we multiply 59, by A, and 4é, 
by Ax Then 


int { (ts cos $+A; cos $—)z sin 
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The quantity in brackets must vanish identically if G # 0, hence 
tan = (Ai +t — (4.3) 


This holds for any form of the function G(t), and has been obtained by 
Fried (2). 

We can solve Eq. 4.3 for sin ¢ and cos ¢, and substitute into the 
integrals for ~ and #7. The result is two transcendental equations to 
be solved for \; and dy. If we assume constant thrust, and let 

= (BA, +1-— Mo)? + XP? = PAP + 


we get 
= — (BA2/a) In + pra? + 
— (aX, + + — 
+ gte)/v.] = In [(Xo + 1 + — ws)/(Xs + BM) ] > (4.5) 
— — we)/a} in + 
+ — (aX, + + — 


With £ = & and » = 0, these are the equations to be solved simul- 
taneously for 4, and Ay. The resulting value of », is then given by 


(B/v.) (ne + 4gte®) = Xo — Xo — (gs + gts) (4.6) 


For a numerical example, we used a hypothetical rocket burning 
94 per cent of its original weight, with 


v, = 104 ft./sec., B = 0.0064 sec.-', ¢, = 146.875 sec. 
g = 32.0 ft./sec.2, pw, = 0.06, &, = 25,856 ft./sec. 


We found fA; = — 0.1420, Br. = 0.5615. The maximum value of m 
came out to be 43.64 miles. 

Remembering now the conclusions of part a, we ask what would 
happen if we do not require 7, to be zero, but instead allow it to be free, 
and maximize 9. + (*/2g) = H, while holding é, fixed. Forming 
5H + dé, = 0, and proceeding as before, we get 


tan @ = (Hj + gts — gt)/d. (4.7) 
This is the same as Eq, 4.3, if we let 
Ar = (mo/g), = (A/g). 


Equations 4.5 still hold for constant thrust except that we do not 
require 7,» = 0, but instead require that 4, = Ag. 
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We gain considerably from not requiring 7, = 0. Using the same 
numerical values as before, we find 
BA = 3.586, wo = 6265 ft./sec., 
n» = 34.94 miles, H = 151.1 miles. 
The gain in maximum altitude is 107.5 miles, a distance in which a 


freely falling body would gain slightly over 6000 ft./sec. in vertical 


velocity. 
The values of tan ¢@ and ¢ are 


tan @ = (342.7 — #)/560.3, ¢ = — 1/560.3(1 + tan? ¢). 


Initially, @¢ = 31.50°, ¢ = — 0.0744°/sec.; at the end of burning, 
@ = 19.30°, 6 = — 0.0910°/sec. 


c. Modifications of the Restricted Problem Produced by a Round Earth with 
an Atmosphere 


If we are dealing with a round Earth with no atmosphere, we begin 
by modifying the equations of motion during burning. For purposes of 
determining the optimum, rather than any specific, trajectory, we can 
surely assume that the motion lies in a plane, and modify the equations 


of motion to 
r—d(r°6)/dt = Gcos¢, *# — = Gsing — (gp*/r*). (4.8) 


In these, we are using circular coordinates in the plane, with r initially 
equal to p, and are using g as the acceleration due to gravity at the 
initial value of 7. @ is the angle made by the rocket axis with the local 
horizontal. 

We assume that the rocket has a first burning period, at the end of 
which it is established in an intermediate orbit. When the intermediate 
orbit is tangent to the final orbit, a second burning period, which we 
shall assume is instantaneous, occurs, at the end of which the rocket 
is in its final orbit. We assume that the total fue! consumption in the 
two burning periods is fixed, but that the partition between the two 
stages is unknown. We wish to determine this partition, and the 
function ¢(#), required to put the rocket into a circular orbit of radius R, 
with a maximum value of R. 

First, we wish an expression for the fuel consumed during the second 
burning period, as a function of the final orbit and of conditions at the 
end of the first burning period. If a subscript b denotes a value at the 
end of first burning, we have* 

8 See any discussion of central force orbits, such as Page, ‘Introduction to Theoretical 
Physics,’ New York, D. Van Nostrand, 1935, Chap. I. 


a 


Sept., 1958.] OptimMuM TRAJECTORY OF A ROCKET 183 


h = = 1476s, (4.9) 
— U = + — = 37 + — (g0*/r), 


along the intermediate orbit. If the final orbit is to be a circle of radius 
R, the value of ¢ in the intermediate orbit must be zero when r = R. 
One readily verifies that 


R = [ge* + — (4.10) 


which connects the desired R with conditions at the end of first burning ; 
one also finds that the velocity in the intermediate orbit when r = R 
ish/R. The required velocity in the final orbit is (gp?/R)*/*, so that an 
increment of velocity AV is required : 


AV = (gp?/R)*? — (h/R). 


Let yp, be the ratio of the mass at the end of first burning to the initial 
mass, and let yu, be the ratio of the mass at the end of second burning to 
the initial mass. Then 


AV = (gp*/R)'? — (h/R) = In (uo/)). (4.11) 


Equation 4.11 is the desired relation involving fuel consumption during 
second burning. 

The most important difference between Eqs. 4.8 and the equations 
in rectangular coordinates is probably the term ré*, which becomes 
slightly larger than g. However, it does so only near the end of burn- 
ing, when G is probably large. The next most important difference is 
probably that between g and (gp*/r*). Hence, a good zero-order 
approximation (fo, 40) to (7, 6) should be the solution to the system 


= Gcos¢, = Gsing — g. (4.12) 
Hence, we can use Eqs. 4.2, if we replace § by pé) — pw, and n by ro — p, 
using w as the initial angular velocity produced by the Earth’s rotation. 
We assume that G(¢t) is known, and that G is reduced to zero at t = th, 


the end of first burning. We vary ¢(¢) and ¢,, but as ¢, is varied, we 
do not vary G, as we did in the burning-time problem. Hence, 


te ) 
= f — t)G cos + Fordte, 


= G cos + Fordto, (4.13) 
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Writing Eq. 4.10 as 


R = to, 6s), 
and 
= + + Ri,56s, 


we have, upon setting 6R = 0 


{R.(t, — t) cos + Ry cos @ — sin 
+ [Ratos + Rat oo + jot, = 0. (4.15) 
Since 6R = 0, we have, upon taking the variation of Eq. 4.11, 


— cos — (ri#/p) sin $}G 
+ + — = 0, (4.16) 


in which , is the burning rate of fuel divided by the initial mass, evalu- 
ated at time ft). 

Multiplying Eq. 4.16 by 4, adding to Eq. 4.15, and equating the 
resulting coefficients of 6¢ and 6t, to zero, we have, for any form of G(t) : 


tan @ = [R, + oi(t, — 


4.17 


with 
= Ry + = (Rin/p) + Ro). 


The form of tan ¢ in Eq. 4.3 still holds, with proper choices for ),, 2. 

Equations 4.5 and 4.6 still hold for values at the end of burning, 
assuming constant thrust. We now have five unknown constants, 
namely 7», 7», 6, 4, t, which must be chosen to fit the five equations, 
namely Eqs. 4.5, 4.6, 4.11, and the second of Eqs. 4.17. The solution 
of these equations is undoubtedly simpler than determining an optimum 
trajectory purely by trial and error, but is still a difficult computing 
task. Hence, we do not give a numerical example. 

The added complexity of the orbit problem in the presence of a 
round Earth does not come primarily from the equations of motion 
during burning, but rather from those considerations of the motion 
after burning which cause the first burning time to be unknown. 

It is interesting to note that the form of ¢(¢) is unaltered, to good 
approximation, by the use of a round Earth instead of a flat Earth, 
provided that we interpret ¢ to be the angle with the local horizontal. 
This is similar to the conclusion reached by Fried and Richardson (1) 
for the maximum range problem. 
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If one wishes to include them, the effects of the atmosphere and of 
the terms neglected in Eqs. 4.12 (centrifugal acceleration, variation of 
gravity, etc.) can all be considered together in a perturbation process 
like that of Sec. 3.c but more complicated. However, the reasons 
which led to results which are quite insensitive to the trajectory when 
we are near the optimum are still valid, and we can use a simple 
process with quite adequate results. The exact process used mighr 
depend somewhat upon circumstances, but a typical one might beast 
follows : 

We can start with the partition of fuel between the two burning 
times which is given by considering instantaneous burning, and assume 
that the consumption during the second burning period is fixed. We 
also compute the intermediate trajectory, and take as the conditions 
at the end of first burning that the rocket is to be established in the 
intermediate trajectory, with one of the variables, say 7,, being maxi- 
mized. Allowing some extra fuel to compensate for losses to gravity 
and drag, and using the shortest burning time allowed by engineering 
considerations (which we assume is far above the optimum; if not, we 
might fix r, — pat about 10° ft. as a good approximation to the optimum 
rs), we calculate a vacuum trajectory using Eqs. 4.5 and 4.6 for values 
at the end of burning, determining the parameters \, and }, to satisfy 
conditions at the end of burning. Thus, we have only two unknowns 
instead of five. 

Finally, to allow for aerodynamic effects, we carry out a perturba- 
tion process similar to that of Sec. 3.c. Instead of using rectangu- 
lar coordinates, we would use curvilinear coordinates, s and u, say, 
with s being arc length along the vacuum trajectory and u being the 
perpendicular distance from s. Use of this coordinate system intro- 
duces acceleration terms other than 4, i, but we already have such terms 
anyway because of the Earth’s curvature. All such acceleration terms 
are small compared with G, and can be put into the perturbation. 
Instead of maximizing the range, we would maximize & subject to 
= 0. 

We expect the effects of the atmosphere to produce results like those 
obtained earlier, namely, causing a rotation of the axis, during the first 
few thousand feet of travel, at a rate required to keep the angle of 
attack small. 


5. SUMMARY 


In Sec. 2 of this paper, we examined the nature of the mathematical 
problem posed by optimizing a rocket trajectory, paying most attention 
to the possibility of auxiliary or ‘‘iso-perimetric’’ conditions. We con- 
cluded that no more than four auxiliary conditions in two-dimensional, 
or six in three-dimensional, motion can be specified, in addition to 
giving some function of the trajectory a stationary value. 
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In Sec. 3, we studied rather thoroughly the problem of maximizing 
the range, including aerodynamic effects. Contrary to a popular belief, 
the rocket should not have low acceleration in dense air and then increase 
its acceleration as it reaches high altitudes. Instead, it should have 
high acceleration in dense air, and should reach the end of burning at 
a comparatively low altitude, such that the loss of range produced by 
drag after burning is a few per cent. The altitude at burn-out should 
be, typically, in the range 50,000 to 100,000 ft. As a rule of thumb, 
we find that the drag should equal the weight under optimum conditions. 

Also contrary to a popular belief, the rocket does not climb almost 
vertically until it reaches thin air, so as to reduce drag. According to 
the theory of this paper, it is true that it starts out rather steeply, and 
then rotates until it is pointed in the desired direction. However, the 
rotation to (almost) 45° is accomplished in a relatively short time. 
By starting out steeply, the rocket gains a little altitude, but rotates 
within a short time in order not to have too much vertical momentum, 
and at a rate which in its details is governed by the need of keeping the 
induced drag small. 

Because the distance travelled during burning is so short for the 
optimum trajectory, effects of the Earth’s curvature during burning are 
small. We may, however, have to consider effects of curvature after 
burning in determining the correct boundary conditions to apply at 
the end of burning. 

In a vacuum, the tilt of the rocket axis should be constant in order 
to achieve maximum range. In air, as indicated above, there is a 
relatively short period during which the axis rotates from almost vertical 
to the desired direction. 

In Sec. 4, we studied the problem of putting a rocket into an orbit 
most efficiently. Since the difference between the air and vacuum 
trajectories should be about the same for this problem as for the maxi- 
mum range problem, we studied only the problem in vacuum. How- 
ever, we outlined a procedure for determining the proper correc- 
tion in air. 

In a vacuum, the tangent of the tilt angle should be a linear function 
of time, so that the rocket should rotate very gradually toward the 
horizontal. This holds for any time variation of the thrust and mass, 
holds exactly for a flat Earth, and holds to good accuracy if we include 
the Earth’s curvature. It is important that the rocket coast during 
most of the flight before it reaches its orbit. 

From the standpoint of applications, an important conclusion is 
that the efficiency is very insensitive to the trajectory. In one example, 
we calculated the optimum trajectory in vacuum, then in air, and 
finally calculated the range which would actually result if the rocket 
followed the vacuum trajectory. The error caused by using the vacuum 
trajectory was about 2 per cent. In another example, we found that 
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the loss caused by using a constant rather than the optimum thrust 
function was only 2 per cent. 

We conclude that there is little need for a more elaborate theory 
of the optimum trajectory than that which we have presented, although 
the details of any specific application may differ considerably. 
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New Compounds Protect Outdoor 
Storage of Bulk Materials.—A new 
compound that blankets and protects 
outdoor storage piles of bulk materials 
against loss from wind or rain erosion 
has been developed by Johnson- 
March Corp., Philadelphia, Pa. The 
spray-applied solution contains speci- 
ally developed additives for maximum 
coverage and penetration. 

The solution, called Compound SP, 
is a blend of synthetic, organic, long 
chain polymers in a water base that 
form a thin crust which is tough, 
durable and highly resistant to 
weathering. In a series of heavier 
concentrations adaptable to a wide 
range of material composition, parti- 
cle size and climatic conditions, they 
are suitable for coating stock piles of 
Bauxite, carbon, chromite concentra- 
tions, coal, Fluorspar, ilmenite, silicon 
carbide, sulphur and many other 
types of materials in outdoor storage 
piles. 

The formed crusts are said to reduce 
the serious economic losses suffered 
each year through the disappearance 
of large tonnages of valuable material. 
In addition, community or plant dust 
nuisances are abated. Protection 
lasts for periods in excess of a year. 

Compound SP is described as being 
inert, unreactive and non-toxic and 
without adverse effect on the burning 
qualities of coal or the subsequent 
processing of sulphur, metal ores, 
carbon, etc. 

Application of the compounds is 
said to be easily accomplished with 
any type of spraying equipment from 
the small portable type up to com- 
mercial orchard spraying equipment. 
The special additives are included to 
provide proper wetting, penetration 
and coverage of the surfaces of even 
the smallest particles on the storage 
pile. 

The solution is supplied in a number 
of different grades ready for spraying 
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at the rate of one gallon of spray per 
100 feet of surface. When desired, a 
color indicator is included at no extra 
cost in Compound SP to show where 
treatment has been made. 


New Ammonium Perchlorate Plant. 
—HEF, Inc., has launched construc- 
tion of a multi-million pound per year 
plant for production of ammonium 
perchlorate, an oxidizer used in solid 
propellants for rockets and missiles. 
The new facility is adjacent to the 
Columbus, Mississippi, plant of the 
Hooker Chemical Corporation of 
Niagara Falls, N. Y.—one of the two 
parent firms. Foote Mineral Com. 
pany of Philadelphia is the other 
owner of the jointly held subsidiary. 

Formation of HEF, Inc., occurred 
earlier this year. After a preliminary 
study of the propellant market, 
Hooker and Foote formalized the'r 
intent to work together in the solid 
fuel propellants field by announcing 
on February 20th that the new 
subsidiary had been formed. The 
HEF, Inc. plant will utilize a process 
developed jointly by the research 
departments of Foote Mineral and 
Hooker Chemical. 

Perchlorates have figured promin- 
ently in recent discussions involving 
rocket and missile propellants. 
Foote has researched certain areas of 
the propellant field, whereas Hooker 
is already a manufacturer of sodium 
perchlorate and potassium perchlor- 
ate. The joint venture is a bid for 
the larger market which is expected 
to develop as a result of the present 
interest in solid propellant systems. 

Hooker Chemical Corporation is 
the largest manufacturer of sodium 
chlorate, the basic raw material for 
perchlorates. Presently Hooker is 
planning to spend one million dollars 
for its third major expansion of 
sodium chlorate production capacity 
at Columbus. 
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A NOTE ON THE OPERATIONAL SOLUTION OF LINEAR 
DIFFERENCE EQUATIONS * 


BY 
ELIAHU I. JURY AND FRANCIS J. MULLIN 


SUMMARY 

The z-transform is applied to the solution of linear difference equations with con- 
stant coefficients; boundary conditions are considered in the resulting solutions and 
certain theorems related to the z-transform are applied to obtain qualitative informa- 
tion about the solution. The z-transform is used to solve simultaneous difference 
equations with constant coefficients and the modified z-transform is used to solve 
difference equations whose coefficients are periodic functions. Finally, the 2z-trans- 
form is applied to the solution of difference equations with periodic coefficients. 

The methods of solution are illustrated with various examples which indicate the 
straightforward procedure for applying the z-transform and the modified z-transform 
to the solution of difference equations. 


INTRODUCTION 


The z-transform (1,2)? and its modified form (3) have been success- 
fully applied to the analysis and synthesis of sampled-data control 
systems. The sampled response of such systems can be described by a 
difference equation with constant coefficients (z-transform) while the 
continuous response can be described by a difference equation whose 
coefficients are periodic functions (modified z-transform). In view of 
this, it is evident that the z-transform and modified z-transform should 
be useful mathematical techniques for the solution of difference equa- 
tions. It is the purpose of this paper to investigate this in detail. 

Transform techniques for the solution of difference equations are 
not uncommon. Gardner and Barnes(4) have used the jump function 
and the Laplace transform, Fort (5) has used the Dirichlet Series trans- 
form, while Neufeld (6) has solved linear mixed difference-differential 
equations with a transform method. Stone (7) has used the General- 
ized Laplace transform to solve difference equations and Brown (8) 
has discussed the theoretical considerations of the method suggested 
by Tustin (9) for approximating differential equations by difference 
equations. Thus, the use of operators and transform methods for the 
solution of difference equations is not new and the z-transform can be 
considered as an extension to these methods. 

* This research was supported by the United States Air Force through the Air Force 
Office of Scientific Research of the Air Research and Development Command, under Contract 
No. AF 18(600)-1521. Reproduction in whole or in part is permitted for any purpose of the 


United States Government. 

1 Electronics Research Laboratories, Department of Electrical Engineering, University 
of California, Berkeley, Calif. 

2 The boldface numbers in parentheses refer to the references appended to this paper. 
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SOLUTION OF LINEAR DIFFERENCE EQUATIONS WITH CONSTANT COEFFICIENTS 
Consider the following difference equation : 


a; and b; are constants, Y, is the forcing function and is known for all 
values of n. It is desired to find X,. X, and Y, are both functions of 
the parameter ¢ but are defined only for = m where 7 is an integer equal 
to or greater than zero.* Let X, and Y, be represented by an impulse 
(delta function) occurring at ¢ = m. With this representation for X,, 
and Y, one may take the z-transform (1, 2) of Eq. 1 and solve for X (z) ; 
X, can then be found by performing the inverse z-transformation. For 
the purpose of completeness it should be noted that we are defining : 


a[ F(n)] F(s) = F(n)s-" (2) 


k-1 
aLF(n + k)] = = Fee) -> (3) 
n=0 i=0 
With these definitions the z-transform of Eq. 1 is 
i—1 


Solving for X (z) one obtains 


X(z)= i=0 i=0 j=0 i=0 j=0 


(5) 


i=0 
The inverse z-transformation can be accomplished using the Cauchy 
Contour Integral (see Eq. 16), the power series expansion method in 
which X (z) is expanded in powers of z~', or with the aid of z-transform 
tables (2,11). It is seen from Eq. 5 that the first p values of X, can 
be arbitrarily chosen to meet given initial conditions; if one neglects 
the initial conditions on both X,, and Y,, a particular solution is obtained. 


Example 1 
Find X,, which is a solution of the difference equation 


— + 0.632 X, = 0.368 + 0.264 Y,,. (6) 
3 In this paper it is assumed that all quantities are 0 for n < 0. 
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The z-transform of this equation when solved for X (z) gives 


(0.368 z + 0.264) Y(z) 4 2Xo + 2(X, — Xo — 0.368 Vo) (7) 
2? — + 0.632 2 — + 0.632 


X(z) = 


If Y, = aconstant = k for all nm > 0, then Y(z) = kz/(z — 1). Using 
a contour integral for inversion we find: 


X, = k[1 — (0.632)*/? cos (0.89 m) — 0.213 (0.632)"/? sin (0.89 m) 
+ (0.632)"/*(1.282 sin (0.89 (m + 1)) 
+ 1.612 (X, — Xo — 0.368 k) sin (0.89 m)]. (8) 


The power series expansion method gives 


X(s) = Xo + + — 0.632 Xo + 0.632 
+ (0.368 X; — 0.632 + 1.265 k)z-* +--+. (9) 


In this equation the coefficient of z-* equals X,, form = k. Needless to 
say, Eq. 9 is identical to Eq. 8 for the proper choice of n. If all initial 
conditions are neglected, then the two forms of the particular solution are 
X, = — (0.632)*/2 cos (0.89 nm) — 0.213 (0.632)*/2 sin (0.89 (10) 
and 


X,(z) = 2° + 0.368 2 + 1.02% + 1.402% (11) 


The particular solution is shown in Fig. 1. 


Fic. 1. Particular solution of Eq. 6 with Y, = a constant k. 


Example 2 
Find X,, which is a solution of the difference equation 


144 
Xn 86 80 895 
4 
2 3 5 6 7 8 
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Following the above procedure we obtain 
3 1 
(g — 

10 X (10 X, 5 (10 5X. + 2 Xo)z 


X(z) = 


102 —5e+22+3 (13) 
and inversion using the power series expansion method gives 
X(z) = Xo + + 0.1(— 3 Xo 
+ 0.111 — —4X1 4+ (14) 
The particular solution for this equation is 
X =0.1 2-*+0.25 2-§+0.405 2-7 +0.543 2-1+0.665 2-8+----. (15) 


This solution is shown in Fig. 2. 


Fic. 2. Particular solution to Eq. 12. 


The contour integral for X, is 


X,= J X (z)2"—"dz. (16) 


It is interesting to note that if X, is to be evaluated using this integral, 
it is necessary to know the roots of the denominator while this is not 
required if one uses the power series expansion method. Also, when 
using the power series expansion method, two theorems which have been 
discussed by one of the authors (2) can be advantageously applied. 
These are 


(a) Initial Value Theorem 


z-—1 


Xo = lim X(z) if X, = 0 for all nm < 0. 


(6) Final Value Theorem 


X, = lim a X(z) if the final value exists. 


Xn 
665 
543 
25 
J 
2790 
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The stability requirements for a.bounded solution are simply that 
the roots of the denominator lie within the unit circle in the z-plane. 
Single roots on the unit circle indicate conditional stability and double 
roots mean instability. It is to be understood that the final value 
theorem as given above applies only if a finite final value exists and 
does not constitute a stability test in itself. The test for the location 
of the zeros of the denominator can be performed using the Schur-Cohn 


criterion (12) or, if one uses the bilinear transformation z = ramey then 
the Hurwitz criterion (12) can be applied to the numerator of the result- 
ing polynomial in w. 


' SOLUTION OF SIMULTANEOUS LINEAR DIFFERENCE EQUATIONS WITH 
CONSTANT COEFFICIENTS 


Consider the following set of simultaneous linear difference equations 


i=0 i=0 i=0 
i=0 i=0 


It is desired to find X, and Y, when F, and G, are known for all values 
of n. To find X, and Y, we take the z-transform of the two equations, 
solve one of the transformed equations for either X(z) or Y(z) in terms 
of the other, and substitute this into the remaining transformed equa- 
tion. This equation can then be solved for one of the unknown trans- 
forms. With one of the unknown transforms determined the second 
can be obtained from either of the original transformed equations. 
Following this procedure for Eqs. 17 and 18 gives 


i=0 j=0 i=0 7=0 

= - = Fj) (19) 


- - 0) 


he 
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From Eq. 19 we obtain 


r 


? 
a;z* 


i=0 


and substituting this into Eq. 20 gives an equation containing Y(z) as 
its only unknown. From this equation we can find Y(z) and knowing 
Y(z) we can obtain X(z) from Eq. 19 by substituting the correct expres- 
sion for Y(z). The inverse transformation will give X, and Y,,. 


Example 3 


Determine the solution of the following set of difference equations 
+ — Xa + Vari — = 2.0 (21) 

+ 3X, + — = 5.0. (22) 

Taking the z-transform of Eq. 21 and solving for X(z) gives 


1 
2+5z—1\z 


X(z)= +2°Xo+2(X1+5X0+ Yo) — (z—3) v(e)). (23) 


The z-transform of Eq. 22 is 
(2+3)X —2Xo+ (42*—1) Y(z) (24) 


Substituting Eq. 23 into Eq. 24 and solving for Y(z) gives 


522+232—11 
z—1 


4 4 Vi) — 3X0~ Xi) 3X1) (95) 
202° — 62?— 52+10 
With the aid of Eq. 25 we can rewrite Eq. 23 as 
ag 
1 42*+20s*— 62*—52+10 


4z4+202*— 6z*—52+-10 


Y()= 


X(z)= 


Equations 25 and 26 represent the transformed solutions of the given 
difference equations and inversion will give the desired values of X, 
and Y,. On carrying out the inversion for the particular solutions 
we obtain 


X,(z) (27) 
Y,(z) =Oz° + +--+. (28) 


4X +23(20X9+4X;) +27( 2Xo- 12 Yo+4 Y;) +2( 2Xo- Xi- Yot 12 Y;) (26) 
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Needless to say, this method is not restricted to two simultaneous equa- 
tions; it can be applied to any number of equations by the obvious ex- 
tension of the above procedure. 


SOLUTION OF DIFFERENCE EQUATIONS WHOSE COEFFICIENTS ARE 
PERIODIC FUNCTIONS 


Consider the following type of difference equation 


q 
(n+p) (n+p +1) i=0 i=0 
fi(t) and g,(t) are arbitrary periodic functions of the independent vari- 
able ¢ with periods equal to 1; X, = X(#)|.-. and Y, isa known function 
of t which is defined only for ¢ = m where a is an integer equal to or 
greater than 0. Again, let X, and Y, be represented by impulses oc- 
curring att = nm. Because of the periodicity of f;(#) and g,(t), we may 
replace f;(¢) and g(t) by f:(m) and g;(m) where m is a parameter which 
varies from 0 to 1. With this change of variable the actual time scale 
becomes t = (nm — 1 + m) and Eq. 29 becomes 


X(n+p+m) = flm)X + gsm) Vay (30) 


Let us define the modified z-transform of X(m — 1 + m) as 


— 1+ m)] = X(z, m) 
and 
+ k +m) ] m). 


With these definitions one can take the modified z-transform of Eq. 30 
and obtain 


Both X(z) and X(z, m) are unknown in this equation. However, X(z) 
can be determined by noting that X, = X(n — 1+ m)|ma:. Thus, 
substituting m = 1 into Eq. 30 gives a difference equation with constant 
coefficients and X(z) can be determined using this equation. The 
resulting expression for X(z) can be substituted into Eq. 31 which can 
then be solved for X(z,m). This would finally give 


? 
fi(m)z* 


? 
1 


I. 

| 
LS 
| 
Y(z) a 
‘= 
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X(n — 1 +m) can then be found using the contour integral since m 
is a constant during the integration process, or with the aid of modified 
z-transform tables (3, 11) or by using the power series expansion method. 
If the power series expansion method were used (11), one would obtain 


2X (z,m) = po(m)z° + pi(m)z-' + + ps(m)z* + ---. (33) 


The coefficient of z-* represents X(t) for k < t < (k + 1) when m is 
varied from 0 to 1. Again, it is apparent that no knowledge of the 
roots of the characteristic equation is required to obtain a solution if 
one uses this method for inversion. The stability requirement for this 
type of equation is the same as that for a difference equation with 
constant coefficients. However, the initial value and final value 
theorems become: 


(a) Initial Value: 


m—0 


s—1 
Final Value: 


lim X(n —1+m) = lim : - ly (z,m) if the final value exists. 
s—1 


O<m<i 
If the final value is independent of m, then there is no steady state 
“ripple’’; if it is not independent of m, then the steady state “ripple’”’ 
can be determined by allowing m to vary from 0 to 1. 
Example 4 
Find X(t) which is a solution of the following equation 


X(t) Xay1(2.582 — 1.582 — #) — X,(1 — 
(m+1)<t< (m +2) + + — 1) + 0.418 Y,(1 — e-*) (34) 


assuming that the coefficients of X;, and Y; are periodic. Replacing the 
t in these periodic functions by m gives the following equation 


X(n +1 +m) = Xo4,(2.582 — 1.582 e-™ — m) + X,(1 — e-*) 
+ + m — 1) + 0.418 Y,(1 — e-*). (35) 


The modified z-transform of Eq. 35 is 


(s,m) = (2.582 — 1.582 e-™ — m)zX(z) — (1 — e-*)X(z) 
+ (e-™ + m — 1)zY(z) + 0.418(1 — e-”)Y(z). (36) 


—<— 
H 
\ 
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To find X(z), let m = 1in Eq. 35 and take the s-transform of the result- 
ing equation. Thus, 


= 0.632 + 0.368 +- 0.264 | 


(37) 


or 


+ — Xo — 0.368 Yo) 
2 — z + 0.632 


(0.368 2 + 0.264) Y(z) 


2? — z + 0.632 38) 


X(z) = 


Substituting this expression for X(z) into Eq. 36 gives 


+m—1)+2(—2.0 — 1.368 m-+2.368) + (e-"+06.368 m—0.737) 
s(@—2+0.632) 

42% 9(2.582—1.582 e~™— m)+2X o(—3.582+2.582 e~-™—m)-+ (2.582—1.582 e~™—m) (X :—0.368 Xo) 

2(2?—z+0.632) 


2(e™ 


Y(z) 


X(z, m) = 


X0—0.368 Yo) 
2(2?—2z+0.632) 


+ - (39) 


The particular solution for X(t) is found using the power series expan- 


sion method with ¥(z) = —*—. This is written as 


(e-™ + m — 1)2° + (0.368 + 0.632 
+ (1.632 — 0.632 + 


2X,(z,m) = 


(40) 


where the coefficient of z-* represents X,(¢) fork <t <k+1. The 
final value of X,(¢) is obtained by applying the final value theorem. 
This shows 


Xin = 


X(n — 1 +m) = lim=— 


The particular solution is shown in Fig. 3. 


LINEAR DIFFERENCE EQUATIONS WITH PERIODIC COEFFICIENTS 


We now wish to apply the z-transform to solve difference equations 
with periodic coefficients (10) ; that is, to the following equation 


X + Kai 


(41) 


a, is a periodic coefficient whose period is equal to N; that is, a; = ay4; 
for 1, 2, ---, (N — 1). is a forcing function known for all 
values of m and Y, and X, are to be represented by impulses. To 
solve this equation we shall first form a constant coefficient difference 
equation which describes the solution form = KN, K = 0, 1,2, --+-, and 
then use this result to find the solution for other values of 2. This 


I | 
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Fic. 3. Particular solution of Eq. 34 for Y, = a constant k. 


procedure, in essence, removes the periodicity of the coefficients from 
the analysis. From Eq. 3 we note that 


aL Yew] = Yo+ + 4+ + --- 


= Y,+ + + + - 


With this pattern for X, we proceed to solve Eq. 41. Let » = KN so 


that a, = do. Wecan then write 

X = + Yun 

= O2X = +0201 V xn + 2 V t+ 


= xn + (@n-1 a;) Yrw 
Notice that the last of these equations given above is actually a constant 


coefficient difference equation involving only Xxy and the various 
forcing terms. The z-transform of this equation can be written as 


| 198 U. FL | 
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2% (Xxw(z) Xo) 
= (Gy-1 @o)Xxw(z) + (Gw-1 G1) 
+ (42a) 
or 
+ + (420) 


Knowing Xxw(z) we can find Xxw4:(z) from the first of the above 
equations, Xxwv42(z) from the second of these equations, and so forth. 
The stability condition is simply that |aoaia:---ay1| <1. If 
Y, = 0 for all » > 0, then 


X xw(z) = 


Example 5 
Find X, which is a solution of Eq. 41 with a; = 0.5 forz = 0, 1, 4, 5, 
8,9, --- anda, = — 0.5 forz = 2, 3, 6,7, 10,11, --- and Y, = (— 1)*. 
As indicated by Eq. 420, we first write 


Yen] = 1+2-*+2-§+--- 
— 


Then using Eq. 426 we write 


Xxn(z) = 


0.125 (0.25)(—24) 0.5 zt 


1 
zi—1 zt—1 zi—1 


or with X>, = 0 
— 1.625 


— 0.0625) (s* — 1)" 


Xxwn(z) = 


Inversion using the contour integral gives 


26 
Xue ~~ 


X = 0.5 Xix + Yux 


Since we know 


= 
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it readily follows that 
4K] 


Similarly we find 


1 
- - 
= € 


The solution to this equation is shown in Fig. 4. It should be noted 


Fic. 4. Solution to Example 5. 


that the inverse transformation can be carried out in a straightforward 
manner for this type of difference equation, since the roots of the 
denominator are readily calculated. 


SOLUTION OF SIMULTANEOUS DIFFERENCE EQUATIONS WITH 
PERIODIC COEFFICIENTS 
Consider the two following difference equations where a,, 5,, c, and 
d, are periodic coefficients of period N, and P, and Q, are known for 
allan. (n > 0) 
X P, (43) 


Voss + + (44) 


A solution for these two equations can be obtained following a procedure 
similar to that used for a single equation with a periodic coefficient ; 
that is, the 2N equations which correspond to a complete period of the 
coefficients are written and each equation is substituted into the follow- 
ing equation. This results in two constant coefficient difference equa- 
tions which involve Xxcwn, Yxw and various forcing function terms. 
These two equations can be simultaneously solved for the two unknowns 
using the procedure introduced earlier in this paper. With Xxw and 


X 4x42 2 
X 4K+3 | 
1.455 1.465 1.465 
125 
Xn 
1875 1335 A 
ia. 7 9 13 4 
-05 
~.934 -934 
~1625 -1735 ~1735 
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Yxw known, one can then find and using 
Xky+1 = + + Paw 
Vena: = CoXxw + doVxw + Qen. 


With and determined, one can find Xxyy:and etc. 
The method is straightforward but unfortunately the notation involved 
in solving two equations of period N is too cumbersome to be useful. 
Rather than solve the general problem of period N, a solution will be 
obtained for N = 2 and N = 3, and this should illustrate clearly the 
procedure to be followed for larger NV. 

From Eqs. 43 and 44 with m = NK we can write 


X wx t+boYwxt+Pwx (45) 
Ywxy1=CoX wx (46) 
vx wx) +4: (coX wx (48) 


wx yx t+dicoX wx +dido (49) 


wx +bicoX wx +bido 


If N = 2, Eqs. 47 and 48 become 


= (Gide + + + Vox 
+ + + (47’) 


= (€1@0 + + + dido) Yox 
+ + diQex + Qoxs:. (48’) 


The z-transform of these two equations can be written 


2*[Xex(z) — Xo] = (aide + + + bide) 
+ a:P2x(2z) + + Poxii(z) (51) 


2° Vex(z) — Yo] = (¢:@0 + dico)Xox(z) + (€1b0 + dido) Yox(2) 
ox(z) + + Qex41(z). (52) 


On solving Eqs. 51 and 52 for X2x(z) and Y2x(z) we find 


Yun — Ky) (2*X0 + H2(z)) + + 
— Ki — K,) — KK; 


Vix(2) = (22 — K1)(2*Yo + Ga(z)) + + H,(2)) 
(22 — K,)(#* — Ks) — KK; 


(53) 


(54) 
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(55a) 


Ki = + 


Kz = aibo + Dido (55d) 
Ks = + (55c) 
= @,Pex(2) + + Pexsi(2) (56a) 


(565) 


G2(z) = + diQex(z) + 


Inversion of Eqs. 53 and 54 will give X2x and Y2x and knowing X2x 
and Y2x one can find Xex,: and Yex,; using Eqs. 45 and 46. If N = 3, 
Eqs. 49 and 50 become 


X KsXsx + + (49’) 
= + KsYsx + G;(z). (50’) 


The z-transforms of these two equations, when simultaneously solved 
for X;x(z) and Y3x(z), can be written as 


(z* Ks) H;(z)) + Ke( Yor? G;(z)) 


Xsx(z) = (s* — Ks)(z* — Ks) — KeK; (67) 
_ — Ks)(z*Vo + Gs(z)) + Ki(Xos* + H;(z)) 
Vix(2) = — Ky) KK (88) 


where 
Ks = + + + (59a) 
= + + + (59d) 
Kz = + + + (59c) 
Ks = + Cabido + detibo + dedidy (59d) 


H;(2) = (aed: + b2¢1)Pax(2) + + Psx+2(2) 
+ (debi + b2d1)Qsx(2) + (60a) 


G3(2z) = + does) Pax(2) + + 
+ + d2d1)Qsx(z) + (605) 


Knowing Xsx and Y3x one can find X3x,, and Y3x4; from Eqs. 45 and 
46 and X3x42 and Yix4: from Eqs. 47 and 48. 


one | 
where 

| 
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Example 6 
Find X, and Y, which are solutions of Eqs. 43 and 44 with 


0.5 meven | 0.0 meven 


P, = —Q. = (—$)**. 


For this example N = 2, so we first compute 


d, = (—1)"# 


— 0.5 


Then using Eqs. 55 and 56 we find 
= — 1; Ky. = 0; K; = 73? K,= — 0.25 


0.5 2? 
z? — 0.25 


— 0.5 2 
2 — 0.25 


== 


G2(z) = 


Finally substituting these values into Eqs. 53 and 54 and simplifying 
the resulting expressions gives, with Xo = Yo = 1.0, 


+ 0.25) 
(2? + 1)(s* — 0.25) 


2*(24 + 1.75 2% — 0.375) 
Yox(z) = (2? + 1)(s* + 0.25) (sz? — 0.25) ° 


Xex (z) = 


The inverses of these expressions are then written as 
= + #(-1)* 
Vox = 4(4)** + 20(4)* — #](-1)*. 
Knowing X2x and Y%2x it readily follows that 
= (—1)*(9)** 
Voxs: = ¥o(4)** — — #](—1)*. 
The complete solution for X, and Y, is shown in Figs. 5a and 5d. 


| 

) 

| 
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Fic. 5. Solutions of Example 6. 


SUMMARY AND CONCLUSIONS 


The theory of the z-transform and the modified z-transform has been 
applied to the solution of linear difference equations. These equations 
often arise in systems which employ magnetic amplifiers or digital 
computers as well as sampled-data control systems and the application 
of the transform method discussed in this paper tends to reduce the 
solution of these difference equations to a straight-forward, systematic 
process. The method solution presented is illustrated with a few ex- 
amples to clarify the procedure involved. Further work on difference 
equations is certainly warranted, in particular, in the field of nonlinear 
difference equations about which little is known at the present time. 
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Test Room for Central Air Condi- 
tioning Systems.—A room simulating 
weather extremes and interior heat 
loads for demonstrating central air 
conditioning systems in virtually any 
new or existing multistory building 
to owners, architects, consultants or 
development engineers, has been an- 
nounced by Carrier Corporation. 

Industry representatives who have 
witnessed the versatile new environ- 
mental room report it is the “most 
completely flexible facility of its 
kind.” It is one of a number of 
features built into a new Systems 
Engineering Laboratory at Carrier 
headquarters in Syracuse, N. Y. 

“Occupants of a building still on 
the drawing boards can learn what 
working in their future office will be 
like,’ Charles V. Fenn, Carrier vice 
president, said. “They will actually 


see the air conditioning unit counter- 
act freezing cold and feel it quickly 


remove scorching heat.” 

Temperature of air in the room, its 
velocity, sound level and distribution 
pattern at any point, as well as effects 
of radiation, can also be noted. 
Design theories for new air handling 
apparatus can be proven under typical 
conditions saving hundreds of man- 
hours. It can also be used to check 
architectural plans and interior motifs, 
eliminating expensive construction 
changes after the building has been 
completed, Mr. Fenn said. 

For example, a prominent Chicago 
consulting firm tested various air 
conditioning systems with a novel 
ceiling it was planning for a new sky- 
scraper. Data collected in one day’s 
experimentation provided essential 


CurRRENT Topics 


F. 


information for designing the building 
around the year-round system. 

The room can be made the exact 
size of standard spaces in a commercial 
structure, hotel or hospital. The 
ceiling is adjustable and the outside 
wall can be changed for glass ratios 
occupying 20 to 80 per cent of total 
area, which includes virtually all 
buildings. External conditions rang- 
ing between zero and 120° F. can be 
created to demonstrate downdraft, 
plus gains and losses due to radiation. 
The effect of venetian blinds or 
drapes can also be shown. Air ve- 
locity of better than 100 feet per 
minute has been recorded sweeping 
down the room’s single thickness 
glass wall when outside temperature 
is zero and 75° F. is maintained 
inside. The degree with which differ- 
ent systems are able to neutralize 
this “waterfall effect”’ is easily one of 
the “most dramatic” demonstrations 
made in the room, Mr. Fenn stated. 

For the finale, the air flow pattern 
provided by one of the air conditioning 
systems is vividly outlined by use of 
a fireworks smoke bomb. 

The room permits comparison of 
the six most popular central systems: 
high velocity, induction; high or low 
pressure ceiling diffusion; high or low 
pressure side wall discharge and 
double duct. Radiant panel sections 
are also installed above the hung 
ceiling. 

By means of an intercommunication 
system and large picture window, in 
one wall of the test module, larger 
groups can follow the scientific 
demonstrations. 
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THE CONJUGATE FRAME METHOD AND ITS APPLICATION IN THE 


ELASTIC AND PLASTIC THEORIES OF STRUCTURES 


BY 
S. L. LEE! 


SYNOPSIS 


This paper discusses the extension of the conjugate beam method to the study of 
the deformation of rigid-frame structures. The first part of the paper deals with the 
elastic analysis of such structures while the second part demonstrates the application 
of the proposed method to the deformation analysis of the structures under ultimate 
loads. The relationship between the conjugate frames and the actual structure is 
derived by means of the “generalized’’ theorems of area moment. Numerical ex- 
amples are given to illustrate the procedure. 


NOTATIONS 


A, B, C, Subscripts denoting points on the frame 


EI Flexural rigidity 

H Horizontal reaction component 
h Length 

M Bending moment 

Mp Plastic moment 

m Area of M/EI diagram 


M pr Area of M/EI diagram between points D and F 


P Concentrated load 

r Number of redundant quantities remaining in the structure at 
ultimate load 

s Number of hinges, plastic hinges plus actual hinges 

V Vertical reaction component 

w Uniformly distributed load 


Xpr, Yor Horizontal and vertical distances between points D and F 
2,9 Centroidal distances of m to any point 
Zp, 9p Centroidal distances of mpr to D, defined by (28, c) 


5., by Horizontal and vertical deflections 

6 Rotation of a section or an actual hinge 

o Rotation of a plastic hinge 

<Mr, Summation of moment about point F in the conjugate frame 


for 6,. 
INTRODUCTION 


The conjugate beam method (1, 2)? provides a powerful tool for the 
study of the deformation of beams. As such a tool, it is commonly 


1 Associate Professor of Civil Engineering, The Technological Institute, Northwestern 


University, Evanston, 
* The boldface numbers in parentheses refer to the references appended to this paper. 
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F. 
employed in the elastic range not only to calculate the slope at any 
section and the deflection at any point of a beam under loads but also 
to determine the redundant quantities in statically indeterminate cases. 
The first part of this paper discusses the extension of this method to the 
elastic analysis of rigid-frame structures, the determination of the 
redundant quantities in statistically indeterminate frames and the 
displacement components of any point in such structures. 

The plastic or yield hinge concept furnishes a simple and rational 
approach for calculating the plastic collapse loads (3, 4, 5, 6) or ultimate 
loads of rigid-frame structures and for estimating the deformation of 
the structures under such loads (7, 8, 9). This simple plastic theory is 
based upon an idealized moment-curvature relationship as shown in 
Fig. 1. The second part of this paper deals with the application of the 


CURVATURE 


Fic. 1. Idealized moment-curvature relationship. 


proposed method to estimate the hinge rotations and the displacement 
components of any point in a rigid-frame structure under ultimate loads. 


THE CONJUGATE FRAME 


Consider a continuous segment DEF (Fig. 2d) of the structure loaded 
as shown in Fig. 2a. The relationship between the displacement com- 
ponents 6p, 5p and 6,» of point D and those of point F may be ex- 
pressed by means of the “generalized” theorems of area moment as 
follows 


Op = Or + Mpr 


byp = — OrXdr — Mpr&p (1) 


= b2r — — 


Mds 
mor= J, “ET 
Mxd 
morte = @) 


p EI 


MOMENT 
where 
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and EI denotes the flexural rigidity of the member. In Eq. 1, xpr and 
yor denote the distances between points D and F along the coordinate 
axes. It is obvious from (2) that mp,» is the area of the M/EI diagram 
while Zp and gp are the centroidal distances of mpr with respect to 
point D. The directions of the bending moment M and the displace- 
ment components 8, 6, and 6, are positive as indicated in Fig. 2d. Thus 
clockwise @ and downward 6, are positive while 6, is positive if the point 
is displaced to the right of its original position. 


80 
(b) CON/SUEATE FRAME FOR dy 
FOR OLF 


dyo 
Oo 


CONJUGATE FRAME FOR bq (d) CONTINUOUS SEGMENT 
FOR OEF OF LOROEO STRUCTURE 


Fic. 2. Relationship between conjugate frames and loaded structure. 


Next consider the conjugate frames shown in Fig. 2b and c. Similar 
to the conjugate beam, the conjugate frames are subjected to the loads 
Mpr, plotted on the tension side of the members, and supported at the 
extremities by the shearing forces 6’s and the bending moments 4’s. 
The only exception is that while the directions of mpr, @p and @r in the 
conjugate frame for 6, (D’E’F’ shown in Fig. 26) are parallel to the 
y-axis, the directions of the same forces in the conjugate frame for 
5, (D’E”F” shown in Fig. 2c) are parallel to the x-axis. With the 
positive directions of mp,r, 6’s and 6’s as indicated in the figures, it is 
obvious that the relationship expressed by (1) can be obtained directly 
from the conjugate frames by simply applying the laws of statics. 
Thus (15) can be derived by equating to zero the summation of moment 
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about point D’ in Fig. 2b, and (1c) can be obtained similarly by setting 
the summation of moment about point D” in Fig. 2c equal to zero. 
Evidently (1a) may be derived from either of the two conjugate frames 
by equating the summation of forces to zero. 

The relationship between the elements of the conjugate frames and 
those of the actual structure and the sign conventions mentioned above 
are summarized in Fig. 3a and b for reference. 


ACTURL CON/UGATE FRAMES 
STRUCTURE FOR by FOR by 


QLNOING MOMENT M 
PLOTTED ON TENSION 
OF MEMBERS 


(0) RELATIONSHIP BETWEEN ME mn 


dya 
POSITIVE DIRECTIONS OF POSITIVE OF SHLARING 
OVSPLACEMIENT COMPONENTS 4ORCES ¢ BENOING MOMENTS OW 


TREMITIES OF CONSUGATE FRAMES, 
QISPLACEMENT COMPONENTS 


POSITIVE O(RECTION OF POSITIVE CONCENTRATED LOAOS OW 
WINGE ROTATION CON SUGATE 


(0) H1WGE ROTATION 


Fic. 3. Relationship between conjugate frames and actual structure and 
sign conventions. 


ELASTIC ANALYSIS 


The method of analysis is best illustrated by means of numerical 
examples. Its application in the elastic range is demonstrated in 
Examples 1 to 4. 
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Example 1 


In order to determine the redundant reaction V¢ and the displace- 
ment components of point C, @¢ and 6,¢, in the frame loaded as shown in 
Fig. 4a, the M/EI diagram is first constructed as shown in Fig. 4d. 


la) ACTUAL STRUCTURE (6) CONJUGATE FRAME FOR Jy 


(e) CONJUGATE FRAME FOR by # (4) ounaram 
Fic. 4. Example 1. 


The vertical arrows in the latter indicate the directions of the loads m, 
acting on the conjugate frame for 45, as shown in Fig. 4b while the 
horizontal arrows indicate the directions of the loads m, acting on the 
conjugate frame for 6, as shown in Fig. 4c. Since point A in the actual 
structure is fixed, its displacement components are zero and, conse- 
quently, the corresponding points in the conjugate frames, points A’ 
and A”, are free ends. For the purpose of illustration, both @¢ and 
5.c in the conjugate frames are assumed to be positive although it is 
apparent from the actual structure that @¢ is negative. 

Considering Fig. 36, the summation of moment about point C’ yields 


wa - (5) wm - ($4) =0 


Ve = #P. 


or 


Also, the summation of forces leads to 


which, upon substitution of the value of V ¢ just obtained, yields 


Ph* 


= 
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Next, considering Fig. 4c, the summation of moment about point C” 


leads to 
Pt) — w (4) a0 =0 


dec = 


Example 2 


Consider the symmetrically loaded square frame of constant EJ 
shown in Fig. 5a. The free body of one-half of the structure cut by 


44? hy 


—= 


Me - 

4 
(a) ACTURL STRUCTURE 4b) FREE BODY OF STRUCTURE 
BY AXIS OF SYMMETRY 
Hits 


ye 


2 oD by 


wk 


(d) CONS FRAME 
(a (e) cons 
Fic. 5. Example 2. 


passing a section through the axis of symmetry is shown in Fig. 5) 
and the corresponding M/EI diagram is shown in Fig. 5c. To de- 
termine the redundant quantities H and M¢, only the conjugate frame 
for 6, shown in Fig. 5d needs be considered. It may be concluded from 
symmetry that 6¢ = 0p = zc = 62» = 0 and, consequently, points C” 
and D” are both free ends. Equating to zero the summation of forces 
and the summation of moment about point C”’ yield, respectively, upon 
simplification, 


2Mc + Hh — &Ph =0 
Mc + — }Ph = 0 


the solution of which yields H = 3;P and Mc = Ph. 
To calculate the vertical deflections of points C and D, b,c and 6,p 
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in the conjugate frame for 5, shown in Fig. 5e are first assumed to be 
positive, or downward. Since member AB does not displace vertically, 
the bending moment at any point in member A’B’ should be zero. 
Considering member A’C’, equating the summation of moment about 


3 
point A’ to zero yields 6,¢ = ineEr Similarly, setting the summation 


Ph 
of moment about point B’ equal to zero gives yp = — eanr the nega- 
tive sign indicating that the deflection is upward. 


Example 3 

To determine the reactions on the two-bay bent of constant EJ 
(Fig. 6a), a shorter solution may be obtained by converting the given 
loads into the symmetrical and the anti-symmetrical systems as shown, 
respectively, in Fig. 6b and c. Thus the redundant quantities are H 4” 
and V4” in Fig. 6b and H,’ in Fig. 6c. 

The conjugate frame for 5, corresponding to ADEB in Fig. 6c and 


“/2 
<1 
(a) LOROEO STRUCTURE (4) SVMMETRICKL LORDING 
w/e 4 
w/e 
% og 
oft 


(0) ANT1- SYMMETRICAL LOR DOING (@) CONJUGATE FRAME FOR 
CORRESPONDING TO ROEB (RIE 


THe 
we 
(e) CON/UERTE FRAME OR J, 1) CONJUGATE. FRAME FOR be 


ta 
CORRESPONDING 70 ROE CAIE CORRESPONDING 70 ADE IW CHIE (6) 
Fic. 6. Example 3. 
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the M/EI diagram are shown in Fig. 6d. Taking the summation of 
moment about A’ eliminates both 6,4’ and @s’ and equating it to zero 
yields H,’ = fe The conjugate frames for 6, and 6, corresponding to 
Fig. 6b and the M/EI diagram are shown respectively in Fig. 6e and f. 
The vertical arrows are associated with Fig. 6e, the horizontal ones 
with Fig. 6f. Observe in this case that the displacement components 
of point E are zero. Setting the summations of moment about point 
A” in Fig. 6e and point A’” in Fig. 6f equal to zero leads to, after the 

3 
elimination of the common factor - 
4H,” 4V + = () 
$H,” 4V 4" + = 0 


solving which yields H,” 
Combining the results of Fig. 6b and ¢ yields the required reactions 


13wh 13wh .. 4wh wh 
Va => 


Hs = 28 336 


Hy, Vz= 


CLOSED FRAMES 
Although symmetrically loaded closed frames may be successfully 


a °°": 


treated in the manner demonstrated in Example 2, the following rela- 
tionships are useful in dealing with non-symmetrical cases. 

If a closed frame is cut by a section at any point, the preservation of 
continuity requires that the displacement components of both ends of 
the cut be identical. Applying (1) to the two end points and observing 
that xpr and ypr become zero in this case, the resulting expressions 
take the forms 


Mpr = 0 
Mprip = 0 (3) 
Mpriip = 0. 


Since the section may be taken at any point around the closed frame 
and, furthermore, it is immaterial where the summations of moments 
are taken, (3) may be generalized in the form 


m = 0 
= 0 (4) 
my = 0. 


In other words, the summation of the area of the M/EI diagram in a 


: 
+ i 
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closed frame and its first moments about any pair of coordinate axes 
must be zero.* 


Example 4 


To illustrate the application of (4), consider the non-symmetrical 
square frame loaded as shown in Fig. 7a. Cutting the section at corner 


CLOSED FRAME (0) OURGAAM DUE 70 EXTERNE 
FORCES 
Ma 
% 
q Ve is 


a 
QAECRAMS DUE TO INTERNAL REOUNDANT FORCES 
Fic. 7. Example 4. 


A, the corresponding M/EI diagrams are shown in Fig. 7) andc. The 
application of (4), taking the summations of moments about corner C, 
leads to 


44M 4 — uh + 4H 4h — Ph=0 
— HV sh + — = 0 
Ma 4h + 4H 4h 4Ph = 0. 


The solution of these equations yields M, = aes V, = sae and 


311P 
555 ° 


DEFORMATION UNDER ULTIMATE LOADS 


Since the shear at any section in the conjugate frame is equal in 
magnitude and sign to the rotation at the corresponding section in the 


*It can be further generalized that the first moment of the area of the M/EI diagram 
about any axis must be zero. 
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deformed structure, the rotation of a hinge in the latter takes the form 
of a concentrated load in the former as illustrated in Fig. 3c. This is 
true for an actual hinge where the bending moment is zero as well as 
for a plastic hinge when the bending moment is equal to the plastic 
moment of the section commonly denoted by Mp. 

The conjugate frame method can be employed with considerable 
advantage over existing methods to estimate the deformation of struc- 
tures at ultimate loads. After the value of the ultimate load is de- 
termined by any method in terms of the plastic moment of the members 
of the structure, or vice versa, the deformation analysis can be carried 
out as illustrated in previous examples with only slight modification. 
If s is the number of hinges, plastic hinges plus actual hinges, and r is 
the number of redundant quantities remaining in the structure at 
ultimate load, it will be shown that (s + r — 1) independent equations 
can be derived from the conjugate frames expressing the relationships 
amongst the hinge rotations and the redundant quantities. When the 
ultimate load is reached, at least one of the plastic hinges is just formed 
where the hinge rotation is zero. This leaves as many unknowns as 
there are equations from which the former can be calculated. 

The location of the last hinge is best explained by means of numerical 
examples. 


(a) LOROEO BENT AND COLLAPSE 
My EL 
MRE CONSTANT 


4. 
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Fic. 8. Example 5. 
Example 5 


Consider the fixed-end bent shown in Fig. 8a subjected to a hori- 
zontal load P applied at joint B and a vertical load P applied at C. 
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It can be shown‘ that the value of the ultimate load is P = 3 ae . The 
collapse mechanism, with plastic hinges developed at A, C, D and £, is 
shown in dotted lines and the corresponding M/EI diagram is shown 
in Fig. 8b. The conjugate frames for 6, and 6, are shown, respectively, 
in Fig. 8c and d in which the directions of the hinge rotations 4, ¢c, 
op and ¢= are determined from the collapse mechanism following the 
sign convention indicated in Fig. 3) and c. 

Since s = 4 and r = 0, 5 +r — 1 = 3 indicating that three simul- 
taneous equations can be derived from the conjugate frames expressing 
the relationship amongst the hinge rotations. Considering Fig. 8c, 
equating the summation of moment about E£ to zero yields 


M ph 


264 6EI = 0. (a) 


In Fig. 5d, taking the summation of moment about £ leads to 


Mrh _ 


and the condition 6,2 = yields 


M ph* 
3EI 


M ph? 
6EI 


+ = + —— 


M ph 
6EI 


= 0. (c) 


For the purpose of illustration, let us first assume that the last hinge 


is formed at C. In this case ¢¢ = 0 which, upon substitution in (a) and 


— and = respectively. Substituting 


the value of $4 just obtained in (c) yields ¢z = ae The negative 


(b), leads to ¢4 = 


value of $4 indicates that this assumption is incorrect since the collapse 
mechanism requires that ¢4 be clockwise as shown in the conjugate 
frames. For the second trial, assuming that the last hinge is formed 
at A yields = 0, dc = op and = Since the 
directions of all the rotations are consistent with the collapse mechanism, 
the second assumption leads to the correct solution. 


‘This example is taken from (10) in which the deformation analysis is carried out by 
means of slope-deflection equations. 
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To obtain the deflections of points B and C, the summation of 
2 


moment about point C in Fig. 8c yields i,¢ = ‘- and taking the sum- 
mation of moment about point B in Fig. 8d leads to 6.2 = sola The 


positive signs indicate that point C deflects downward and point B 
moves to the right. 


Example 6 


Consider the two-story bent loaded as shown in Fig. 9a. It can be 


demonstrated that the ultimate load* is P = a the collapse mecha- 


nism being indicated in dotted lines with plastic hinges developed at 


2 


(a) TWO-STORY BENT AWD (0) COWLSUGRTE FRAMES 
COLLAPSE MECHANISM. FOR by FOR RED. 
Mp ££ ARE CONSTANT. 


CONSUGRTE FRAMES (d) CONSUCRTE FRAMES 
COMF 


Fic. 9. Example 6. 


points A, B, C, D, E, Fand G. The corresponding M/EI diagram is 
shown in Fig. 95 while the conjugate frames for CDH F and ABFG are 
shown, respectively, in Fig. 9¢ and d. To reduce the number of figures, 
the conjugate frame for 5, and that for 6, are combined in a single figure 
in each case. The vertical arrows are associated with the latter, the 
horizontal ones with the former. As in the last example, the directions 
of the hinge rotations are taken consistent with the collapse mechanism. 


5 See (6) for the determination of the ultimate load. 
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It should be observed that, in accordance with the sign convention 
shown in Fig. 3, the directions of ¢c, ¢r and the loads on member CF in 
Fig. 9c are opposite to those on member BF in Fig. 9d. 

In this case s = 7 and r = 0, hence s + r — 1 = 6 indicating that 
six independent equations can be derived from the conjugate frames. 
In Fig. 9c, = 0,6 Mc, = Oand = 0 yield, respectively, 


5M ph 


én — dc + Ser 7° 


M ph 


ph 


260 


Considering Fig. 9d, Ms, = 0, “Me, = 0 and equating 6.2 = 6.r 
lead’ to 


M ph 


or = 9 


Mrh 
os + bc — 64 — Fonz = 9 (e) 


M ph 
da — + = 0. 


er that the last hinge is formed at D, ¢p = 0 and (d) and 
_ SMeh 35.M ph 11M ph 31M ph 


oa = uae and ¢3 = a It should be noted that the solution of 
(d) and (e) involves very little numerical work due obviously to the 
fact that each equation involves only two or three unknowns. 

To estimate the lateral deflection of each story, taking the summa- 
tion of moment about point D in Fig. 9e yields 


5M ph? 
3EI 


6 For brevity, 2 Mr, denotes the summation of moment about point F in the conjugate 
frame for 

2M ph* 
= — dah — ET and = — — , both assumed positive or toward the 


right. 


6EI ’ 72EI ’ ’ 36EI 
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Taking moment about point B yields 


5. — 


The negative signs indicate that the deflections are toward the left. 


Example 7 


To illustrate the treatment of frames with inclined members, con- 
sider the two-bay gable frame (Fig. 10a) which is subjected to a hori- 


(0) TWO -8RY GRELL FRAME ANDO COLLAPSE MECHANISM 
CAF 


(4) LVACRAM 


CONJUSRTE FRAME FoR Ty (@) FRAME 


(€) CONJUGATE FRAME FOR dy 


Fic. 10. Example 7. 


zontal load P at joint D and two vertical loads 7P and 5P applied at 
joints E and G, respectively. The plastic moment and the flexural 
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rigidity of the members as well as the collapse mechanism are indicated 
in the figure. It can be shown that the ultimate load is P = a with 


plastic hinges developed at joints D, E and F. Under ultimate loads 
the structure is statically indeterminate to the first degree and, taking 
He as the redundant reaction, the corresponding M/EI diagram is 
shown in Fig. 10d. In this case s= 6, r=1 and s+r—1=6 
indicating that six independent equations can be derived relating the 
rotations of the plastic hinges ¢p, ¢z and ¢,, the rotations of the actual 
hinges 64, 0s and @¢ and the redundant reaction H¢. 

In the conjugate frames shown in Fig. 10c, d and e, the directions of 
op, dz, ¢r and 6,4 should be consistent with the collapse mechanism. 
After these four quantities and H¢ are determined by means of four 
of the six equations, 0s and 9- may be calculated from the remaining 
two equations with arbitrarily assumed directions, the correct directions 
being indicated by the signs of the solution. 

In Fig. 10d, } Mc = 0 yields 

Mp 


He = 1.283 (f) 


Substituting (f) in Fig. 10 results in an M/EIJ diagram expressed in 


terms of of Considering Fig. 10c, }Ms = 0 yields 


EI’ 
h 
In Fig. 10e, >Mr = Oand >My = 0 lead to, respectively, 
M 
M ph 


264 + 1.5¢2 dr 0.904 = ud 0. 


Assuming the last hinge to be formed at D, ¢p = 0 for which the 
Mrh 


solution of (g) and (h) yields ¢¢ = 0.189 ——,, dr = 0.191 —— Fr and 
6, = 0.406 

With the values of 64, dp, ¢z and ¢r known, equating the summation 
of forces to zero in Fig. 10c yields 6, = 0.275 ue and the same equi- 


librium condition in Fig. 10d leads to 6¢ = 0.436 ae With all the 


hinge rotations known, the rotation and deflection of any point in the 
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frame can be computed simply by calculating the shear and moment 
at the corresponding point in the conjugate frames. 


CONCLUSIONS 


It is evident that the conjugate frame method, when applied to beam 
problems, degenerates into the conjugate beam method. As obvious 
is the fact that the procedure employed in the last three examples lends 
itself readily to the deformation analysis of beams under ultimate loads. 

Although no examples are given for frames composed of members 
with varying moment of inertia, the analysis of such cases presents no 
new problem other than the determination of the area and the location 
of the centroid of more complicated M/EI diagrams. In cases where 
these diagrams are higher-degree curves, the integrations defined by 
(2) may be carried out or they may be approximated by summations 
with the curves replaced by step functions. The above comments 
apply equally well to frames with curved members. 

In the plastic design of rigid-frame structures composed of structural 
steel members, the ductility of the materials ensures, in common cases, 
that the hinge rotations required to develop the plastic collapse load 
can be sustained by the plastic hinges. Even then a deformation 
analysis should prove to be both wise and necessary in order to judge 
whether the deflection of the structure or portion thereof due to acci- 
dental overloading is within acceptable limit. In the ultimate strength 
design of reinforced concrete structures (11), the small strain capacity 
of concrete turns the hinge rotations into the major design consideration. 
To this end, among others, the conjugate frame method may prove to 
be a useful tool in addition to the ones heretofore available. 
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E2 — M1 MIXING RATIOS IN 2+ — 2+ — 0+ TRANSITIONS * 
BY 
V. R. POTNIS ' AND C. E. MANDEVILLE * 


When a nuclear energy level decays by emission of electromagnetic 
radiation, depending upon the spin sequence, transitions of any one 
of several different multipole orders may be possible. In the case of 
a 2+ —+2+ —0+ spin arrangement, the first transition, 2+ — 2+, will 
often be magnetic dipole or electric quadrupole in character. Such a 
transition is said to be mixed. The gamma ray emitted second in time 
always results from a pure E2 transition, the spin change ranges from 
2 — 0 to 2 + 0, and so has only the value 2. In the case of the mixed 
transition, the relative intensities, J(E2) and I(M1), are respectively 
proportional to the squares of the matrix elements which describe the 
two types of transition. The quantity 6 is defined as 


§ = + VI(E2)/I(M1). (1) 


In a recent report (1),? Sakai has sought to associate certain values 
of 6 with the 2+ — 2+ — 0+ de-excitation of even-even nuclei in certain 
regions of the atomic table. More specifically, he has attempted to 
develop a form of systematics involving the nuclear shell model which 
relates the sign and magnitude of 6 to regions of neutron number which 
precede the closing of the go/2, 411/2, and 413/2 sub-shells. These “islands” 
of nuclei occur for 38 < N < 48, 64 < N < 74, and 110 < N < 124. 
From his consideration of published data, Sakai has concluded that the 
sign of 6 is negative in the first and third islands and positive in the 
second. 

The value of 6 is often determined from spatial correlation experi- 
ments involving the two emitted gamma rays of the 2+ — 2+ —0* 
transitions. Because of the systematic significance which has been 
attached (1) to the value of 4, previously reported measurements of this 
quantity for some even-even nuclei of the third island have been 
re-evaluated. 

The angular correlation function for two gamma ray quanta emitted 
successively in cascade is 


W(6) = 1 + (cos 0) + (cos 8) (2) 


* Assisted by the joint program of the Office of Naval Research and the U. S. Atomic 
Energy Commission. 

! Permanent address: Gwalior (M.P.), India. 

? Assistant Director. 

* The boldface numbers in parentheses refer to the references appended to this article. 
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where P, (cos @) and P, (cos @) are the Legendre polynomials, and the 
coefficients A; and A, depend on the spins and the multipole orders of 
the emitted radiations. The mixing ratio, 4, is related to the coefficients 


A; and A, by the equations, 


1 
= (0.250 + 0.7326 — 0.0778) 
(3) 
_ 0.3268 


Equations 3 are plotted in Fig. 1. The quantity 6 can also be obtained 


A; 
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Fic. 1. Coefficients A, and A, plotted as a function of 6. Values of A: and A, corresponding 
to 6 negative for Os'** are indicated by the open circles. 


from the anisotropy A defined as 


_ W(180°) — W(90°) 


A 


or in terms of A; and A,, 


(5) 


Equation 5 is of course obtained from Eqs. 2 and 3. Equation 5 can 
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further be reduced to 


_ 3 + 8.7826 + 0.7158 


A= 79.9078 + 9.2868" (6) 


From experimentally determined correlation functions, values of A, and 
A,, and of A, can be obtained. 

The observed values of A, and A, are plotted in Fig. 1. For the 
cases of Os'**, Pt!*, and Pt!%*, the observed values (3-5) of Az and A, 
lie upon or near curves calculated for a positive value of 6. The 
values (2) of Az and A, for Os'** could be interpreted as suggesting 
either a positive or negative value of 6. 

The data for the four elements under consideration are summarized 
in Table I, where the positive value of 6 for Os'**, Pt, and Pt! has 


TABLE I. 
Ag 


A 6 


Nucleus Reference Az 


2 —0.02 +0.01 0.18+0.03 0.076+0.023 

Os'88 3 —0.037+0.002 0.35+0.02 0.14 +0.01 +15.942.4 
ye 4 —0.05 +0.003 0.30+0.015 0.099+0.009 +44.6+18.3 
5 0.075+0.004 0.28+0.02 0.269+0.013 +5.48+1.10 


been calculated from Eq. 6. Because of the large error in A for Os!**; 
and because of the uncertainty of sign, no value of 6 for Os'** has been 
calculated from Eq. 6. 

From Table I, it is clear that except for Os'**, the sign of 6 is defi- 
nitely positive for Os'**, Pt! and Pt!*. Thus, of the seven isotopes in 
Sakai’s third island (1) for which the value of 6 is available, three are 
positive, three are negative, and one, that of Os'**, cannot be determined 
definitely from the presently available data. Thus the conclusion of 
Sakai, that the sign of 5 is negative in this third island, is not fully 
substantiated by the presently available data. Whether any system- 
atic relationship between 6 and the neutron number exists, remains 
uncertain. 


(1) Mitsuo Saxal, Report on Even-Even Nuclei, Institute for Nuclear Study, University of 
Tokyo, Japan. INSJ-6, February 24, 1958. 

(2) T. Lrunpoevist AND I. MARKLAND, Nucl. Phys., Vol. 4, p. 189 (1957). 

(3) V. R. Potnis, V. S. Dupey anp C. E. MANDEVILLE, Phys. Rev., Vol. 102, p. 459 (1956). 

(4) C. E. MANDEVILLE, J. VARMA AND B. L. Sarar, Phys. Rev., Vol. 98, p. 94 (1955). 

(5) V. R. Potnis, Indian Journal of Physics, Vol. XXX, p. 375 (1956). 
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AN ELECTROLYTIC METHOD FOR TRANSIENT MIXING 
MEASUREMENTS * 


BY 
MARVIN P. NORIN! 


In certain fluid systems it is necessary to know the transient tempera- 
ture behavior of a fluid volume in the system when the system is sub- 
jected to a change in operating conditions. A nuclear reactor vessel, 
for example, may contain large plenums or volumes of coolant wherein 
the fluid is at a certain temperature. As a consequence of a change in 
reactor operation such as an emergency shutdown, the temperature of 
the fluid entering the fluid volume can rapidly change to some new 
value. The temperature of the incoming fluid can change by as much 
as 300° F. in about one second. In the engineering design of the 
reactor vessel, it is necessary to know the time variation in temperature 
of the fluid at locations near the vessel structure in order to compute 
the ensuing transient thermal stresses imposed on the vessel structure. 
For reactor systems studies the variation of the temperature of the fluid 
leaving the vessel is required so as to evaluate this effect on other com- 
ponents of the system. 

The time variation in temperature of a fluid volume in a flow system 
is usually due primarily to the mixing of portions of fluid at different 
temperatures. That is, the inertial properties of the fluid are para- 
mount, subject, of course, to the geometry of the containing structure. 


CONDUCTIVITY CELL 


ELECTROLYTE 
——> > VESSEL 
MAIN FLOW 
( woter) 
BRIDGE RECORDER 
CIRCUITS 


Fic. 1. Electrolytic method for transient mixing measurements. 


* Work sponsored by Atomic Power Development Associates, Detroit, Mich. 
1 Senior Research Engineer, Heat Transfer & Fluid Mechanics Branch. 
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The thermal properties of the fluid may be neglected. For this case, 
the temperature changes in transient mixing may be simulated by an 
electrolytic method which has been developed at The Franklin Institute 
Laboratories. 

To illustrate the method, consider the diagram of Fig. 1. What is 
required here is the mixing characteristics of the vessel such that the 


Fic. 2. Conducting a transient mixing test on a flow model of a nuclear reactor vessel. 


time variation in temperature of the fluid leaving the vessel can be 
determined for a fast change in temperature of the fluid entering the 
vessel. 

In the electrolytic analogy, a steady stream of electrolyte is sud- 
denly injected into the main water flow and mixes with it before entering 
the vessel. Thus the concentration of electrolyte in the incoming flow 
is rapidly changed from zero to some pre-determined value. The 
electrical resistance of the incoming fluid is sensed by a flow-type 
conductivity cell. 

Another conductivity cell is located at the vessel outlet and senses 
the resistance of the fluid leaving the vessel. The conductivity cells 
are connected in the active arms of bridge circuits. The bridge circuits 
are connected to a multiple-channel recorder which givés a time record 
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of the electrical resistance or conductance at the cell locations. These 
data can then be converted to values of electrolyte concentration. 

Now if the inertial properties of the fluid are paramount in the mix- 
ing process, then the diffusion properties of the electrolyte may be 
neglected in accounting for the changes in concentration. Thus, the 
changes in concentration are closely analogous to temperature changes. 
The mixing characteristics measured with the electrolytic method can 
then be used to predict transient behavior. 

There is no inherent inertial-type time lag in a conductivity cell 
such as that found ina thermocouple. Any time lag is only that neces- 
sary for ‘“‘washing away” of fluid at the electrodes. This can be 
minimized by good design of the cell so as to allow free passage of the 


fluid past the electrolyte. 


re INPUT TO VESSEL (incoming fluid) 


RESPONSE AT VESSEL OUTLET 


CHANGE IN CONDUCTANCE 
MAXIMUM CHANGE IN CONDUCTANCE 


TIME ——> 
Fic. 3. Sample response curve—transient mixing in a nuclear reactor vessel. 


We have used sodium hydroxide solution as the electrolyte. Our 
experience with the method has shown that large changes in electrolytic 
conductance, which allow ease of measurement, can be introduced with 
very small quantities of electrolyte. This is important since for dilute 
solutions electrolytic conductance is very nearly linear with concentra- 
tion and therefore the conductance measurements may be used directly. 
By proper design of the electrolyte injection circuit, it is possible to 
approach a step change in electrolytic conductance of the incoming 
fluid (see Fig. 3). 

The output of the instrumentation will not, in general, be linear over 
the ranges used. A calibrated resistance circuit can be built into the 
instrumentation so that each set of data can be provided with a set of 
calibrated values. It is necessary to use a-c. power for the conductivity 
cells so as to avoid polarization.. We have used a 2500-cps. power 


supply. 
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An experimental setup for mixing measurements in a flow model of 
a nuclear reactor vessel is shown in Fig. 2. A sample response curve for 
a cell location in the vessel outlet is given in Fig. 3 for the input shown. 
In Fig. 3 the dimensionless ratio of the change in conductance to the 
maximum change in conductance is plotted against time. By convert- 
ing to concentration ratios (not necessary in the linear range of the 
conductance versus concentration curve) one can use the curves of 
Fig. 3 to predict the time-temperature behavior at the vessel outlet in 
response to the change in temperature of the incoming fluid. This 
response can also be used to calculate the response to inputs other than 
step changes. Or, another approach would be to design an electrolyte 
injection system which would produce the desired input and measure 
the response directly. 
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Seminar on A Physical Chemical Study of Calf Thymus Deoxyribo- 
nucleic Acid.*—Dr. VERNE NoRMAN SCHUMAKER of the Department 
of Biochemistry, University of Pennsylvania. Philadelphia, Pennsyl- 
vania, addressed the Seminar of The Biochemical Research Foundation 
on Tuesday, April twenty-second. 


Fresh problems, significant and challenging, are afforded the physi- 
cal biochemist by the nucleic acids. Deoxyribonucleic acid (DNA), 
according to recent concepts, is the bearer of the genetic information 
of the living cell ; it may be considered to be a macromolecular blueprint 
containing the genetic information coded by the sequence of purines 


(1/C.) dc/ds 


0 35 
20 
Sow 
Fic. 1. The sedimentation coefficient distribution of calf thymus DNA. The abscissae 
are the sedimentation coefficients corrected to 20°C. in water. Since the concentration at 
which the ultracentrifuge run was made was 0.001%, correction for concentration dependence 
can be neglected. The ordinates are values of the normalized distribution function, (1/Co) 


(dC/dS). Corrections for the radial shape of the centrifuge cell and for the inhomogeneous 
centrifugal field which would distort the distribution function have been applied. 


45 50 55 60 65 


* The work reported in this communication was performed at the Virus Laboratory, 
University of California, Berkeley, Calif. This work was supported in part by a grant from 
the National Science Foundation. 
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and pyrimidines along a doubly stranded helix. Hence, the study of 
the physical chemical properties of this molecule may be of importance 
as a first step toward ‘‘cracking”’ the genetic code. 

The difficulties inherent in measuring the molecular weight as well 
as many other physical properties of DNA are intimately related to the 
tremendous asymmetry of this biological polymer. Hydrodynamic 
interactions between these extended molecules occur even in fairly 
dilute solutions and are pronounced at the concentrations employed 
when using schlieren optical techniques. Hydrodynamic interactions 


INTRINSIC viScosiTy 


Fic. 2. The relationship between the intrinsic viscosity and the sedimentation coefficient 
for calf thymus DNA. The ordinates are the intrinsic viscosities in units of 100 cm./gm. 
The abscissae are the sedimentation coefficients essentially at infinite dilution, corrected to 
20° in water. 


are not significant if extremely dilute solutions are studied. For ex- 
ample, ten years ago ultracentrifuge measurements of calf thymus 
DNA had always been performed at concentrations sufficiently high to 
resolve the sedimenting boundary with the schlieren optical system (1). 
Recently, when concentrations were reduced 50 fold and the sediment- 
ing boundary observed with an ultraviolet absorption optical system, 
the observed properties of the molecule changed markedly. With the 
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schlieren optical system, a sharp spike with a sedimentation coefficient 
of 10 svedbergs had been observed. At the lower concentration and 
employing UV optics, there was now seen to be a whole spectum of 
molecules with sedimentation coefficients ranging between 5 and 60 
svedbergs (see Fig. 1) (2, 3). 

This extreme polydispersity in sedimentation coefficient exhibited 
by calf thymus DNA could conceivably be caused by a corresponding 
polydispersity in either the molecular weight or the shape of these 


(1/Co) dC/dM 


we 2024 
r ais 
10,000) x10" 
Fic. 3. The molecular weight distribution of calf thymus DNA. The ordinates are the 
normalized distribution function (1/Co) (dC/dM). The abscissae are the molecular weights 
X 10-* multiplied by a factor (10,000/K), which is a function of the shape of the molecule and 
should be close to unity. 


particles in solution. In order to decide between these two extreme 
possibilities, the DNA preparation was divided by preparatory cen- 
trifugation into 10 fractions, each fraction containing molecules of 
different sedimentation characteristics. The sedimentation coefficient 
distribution and the intrinsic viscosity of each fraction were then meas- 
ured, and from this information the relation between the sedimentation 
coefficient and intrinsic viscosity was worked out for calf thymus 
nucleic acid. This relation is shown in Fig. 2; the intrinsic viscosity 
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increases as the sedimentation coefficient increases, which is just what 
would be expected if the polydispersity in sedimentation coefficient 
were primarily due to a polydispersity in molecular weight. A decrease 
in the intrinsic viscosity would be found if the polydispersity were due 
to shape factors alone. 

The molecular weight distribution can be calculated from the data 
of Figs. 1 and 2, assuming that the polydispersity is entirely due to the 
variation in the molecular weight. This has been done and is presented 
in Fig. 3. It can be seen that this preparation of DNA has a range of 
molecular weights between 4 and 25 million, with a weight average 
value of 10 million. 


VISCOSITY OF DNA SOLUTION 
Ordinate: SPECIFIC VISCOSITY 
Abscissa: MINUTES AFTER ENZYME ADDITION 


© Lower Bulb Shear Gradient 43 Sec* 
4 Middle Bulb Shear Gradient 92 Sec™' 
oO Upper Bulb Sheor Gradient 145 Sec-! 
o Extrapolated Shear Gradient O Sec"! 


— i L L i i 
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Fic. 4. The viscosity change during enzymatic degradation of calf thymus DNA as a 
function of time. The three lower curves represent the experimentally determined values of 
the reduced viscosity at three different shear gradients. The upper curve is the extrapolated 
reduced viscosity at zero shear gradient. The ordinate is the specific viscosity, and the abscissa 
is time in minutes. 


Occasionally, knowledge of intimate structural detail may be gained 
by physical-chemical studies. Important architectural information 
about DNA has been obtained by a study of the kinetics of the enzy- 
matic degradation of this molecule. In solution DNA is a rope-like 
structure. If this “rope’’ has but one strand, then at each attack by 
the degrading enzyme, it should break, and the result should be a large 
decrease in the molecular weight and, for highly asymmetric particles, 
the intrinsic viscosity. If the DNA ‘rope’ is composed of two strands, 
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held together by weak crosslinkages, then the molecule will not break 
until the enzyme has cleaved both strands at adjacent positions along 
the rope. Hence, a study of the intrinsic viscosity of a solution of 
DNA undergoing enzymatic degradation should afford considerable 
structural information. This study has been made, and the results are 
are presented in Fig. 4. It is obvious from this figure that there is a 
protracted lag period after the introduction of the enzyme before the 


EFFECT OF HEAT ON DNA 
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Fic. 5. The degradation of DNA caused by heating at 100°C. in NaCl solutions of 
various molarities. These unbuffered solutions have pH values of about 6.5 at room tempera- 
ture. The ordinates are the logarithms of the molecular weights. The abscissae are the times 
of heating at 100°C. 


rate of decrease of the intrinsic viscosity attains its maximum value. 
This lag period is characteristic of a multi-stranded structure, and the 
one strand rope model for DNA may be eliminated at once. Detailed 
mathematical analysis of the shape of the intrinsic viscosity vs. time 
curve of Fig. 4 shows that the data are best described by a two-strand 
rope model for DNA and that both of the strands are continuous along 
the length of the molecule, containing few if any preformed inter- 
ruptions (4). 

Additional evidence for a non-interrupted doubly-stranded struc- 
ture for DNA comes from examination of the effects of heat degradation. 
Heating solutions of DNA at 100° C. at various ionic strengths causes 
slow but progressive breakdown of the macromolecule. (See Fig. 5.) 
However, the intrinsic viscosity decreases 50 fold during the first 2 or 3 
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minutes of heating (5). This suggests that the molecular rigidity, 
which gives a highly extended structure possessing a large intrinsic 
viscosity, is due to multiple hydrogen bonding between the two chains. 
Upon heating, the weak hydrogen bonds break at once, and the intrinsic 
viscosity falls markedly. Changes in the molecular weight, however, 
reflect cleavage of strong “‘backbone’’ bonding—probably the 3’,5’- 
phosphate diester linkages. Some information as to how these back- 
bone bonds may be broken is given in Fig. 6. When DNA is heated 
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Fic. 6. The left-hand figure shows the release of the purines adenine and guanine from 
DNA as a function of the time of heating at 100°C. in a neutral but unbuffered solution. 
Paper chromatography was used for the identification and the estimation of the bases. The 
right-hand figure shows the relative optical density increments obtained after heating DNA 
at 100°C. at various pH values for one hour, then precipitating the DNA with acid and 
examining the supernatant liquid. 


HOURS 


at 100° C. at acid or neutral pH values, a great amount of ultraviolet 
absorbing material is released into solution. By paper chromatography 
this material yields the purines adenine and guanine, an indication of the 
sensitivity of the N-glycosidic bond to heating, particularly at acid pH 
values. The splitting of the N-glycosidic bond frees the reducing end 
of the sugar which may then activate the cleavage of the diester phos- 
phate bond. 

In summary, application of physical chemical techniques has given 
a detailed picture of the molecular weight distribution of a sample of 
calf thymus nucleic acid, which extends from } to 25 million, and yields 
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a weight average value of 10 million. The structure of the DNA 
molecule has also been examined and the results obtained are best ex- 
plained by a two-strand rope model. The integrity of the strands and 
the sensitivity of various chemical bonds within the molecules to heat 
treatment have been included in this study. 
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NucLear RapiaTion DETECTION, by William 
J. Price. 382 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1958. Price, $9.00. 


The era during which users of nuclear 
radiation detectors designed and constructed 
their own equipment has passed on. Gener- 
ally speaking, commercial sources now supply 
most of the required instruments, and, in- 
deed, one is likely to become confused by the 
large variety of detectors which are now 
available. 

This book presents, in a most intelligible 
manner, basic scientific information which 
should enable even a novice to determine 
which type of instrument best fulfills his 
particular requirements. On the other hand, 
the experienced research physicist who is 
already thoroughly conversant with the 
principles of operation and properties of the 
various detection devices covered will welcome 
this extremely well-organized collection of 
basic material which has heretofore been 
difficult to extricate from the very extensive 


and widely-scattered literature on the many 


aspects of the subject. Furthermore, funda- 
mental considerations which are important 
in the interpretation of the measurements 
obtained with the various types of detector 
are discussed; this is a matter of obvious 
importance to all who design or use radiation 
detection equipment, or who are dependent 
upon measurements made by others. 

Chapter 1 deals with the properties of 
nuclear radiation, specifically with respect to 
interactions with matter. A brief survey of 
detection methods, followed by a discussion 
of the statistical nature of detecting systems, 
introduces chapters dealing with ionization 
chambers, Geiger-Muller counters, propor- 
tional counters, scintillation detectors, 
photographic emulsions and other detection 
methods (that is, cloud chambers, crystal 
counters, Cerenkov detectors, chemical do- 
simeters, calorimeters). A separate chapter 
is then devoted to neutron detection methods, 
and finally there is a detailed treatment of 
the principles of operation of some of the 
electronic circuits associated with nuclear 
radiation detectors. 

The book is well-written and illustrated, 
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contains an ample bibliography, and is 
logically organized. It will undoubtedly 
prove of great value in the training of students 
of nuclear science and technology, and would 
serve admirably as a text for a course in 
nucleonics instrumentation. There are nu- 
merous illustrative examples with solutions 
interspersed throughout the text, as well as a 
collection of problems at the conclusion of 
each chapter 
In the opinion of the reviewer (who is 
admittedly biased) the omission of a chapter 
on cosmic-ray detectors, which would have 
been most timely, is unfortunate. There are 
only occasional references to the application 
of some of the techniques to cosmic-ray 
experiments, but no consistent policy is 
evident. Nevertheless, Nuclear Radiation 
Detection is an excellent book which certainly 
serves the purposes for which it was written. 
M. A. PoMERANTZ 
Bartol Research Foundation 


THE GyRoscOopPE—THEORY AND APPLICA- 
TIONS, by James B. Scarborough. 257 
pages, diagrams, 6 X 9in. New York and 
London, Interscience Publishers, 1958. 
Price, $6.50. 


While the dynamics of the gyroscope are 
included in most books on mechanics, we 
also find special treatises, such as Klein and 
Sommerfeld’s ‘Theorie des Kreisels,’’ Gray’s 
“Gyrostatics and Rotational Motion,”’ Gram- 
mel’s Kreisel” and Greenhill’s ‘““Gyro- 
scopic Theory,” all of an advanced nature. 
At the other extreme, there are numerous 
elementary books of a popular nature. In 
the intermediate class we have Crabtree’s 
“Spinning Tops and Gyroscopic Motion” 
and “Mechanics of the Gyroscope,” by 
Deimel. Prof. Scarborough’s book is a 
welcome text adding to this latter group, 
with a good analytical and vector treatment, 
with many practical applications on important 
engineering problems. 

The theory of the gyroscope well illustrates 
the principle of the space rate of change in 
angular momentum; and the reactions set up 
are the moments or torques required to cause 
these changes. In the more complicated 
problems we have couplings between the 
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various coordinates defining the system. 
Such gyroscopic couplings are of importance 
in the understanding of the operation of 
technical applications. For these applica- 
tions, as well as an elementary understanding 
of the basic theory, the analytical approach, 
as used by the author, is needed. Of special 
interest is the well chosen group of applied 
problems, such as the top, ship stabilizers, 
the gyroscopic compass, direction indicators, 
the monorail car, applications to machinery 
as well as some astronomical applications. 

The author’s analytical approach to the 
various problems could be considerably 
improved by introducing a more physical 
discussion of the interactions. Thus, in the 
discussion of the rolling hoop, it seems un- 
necessary to introduce a prescribed curvature 
with radius R and corresponding “‘centrifugal”’ 
force mv* to account for the lateral space-rate 
of change of linear momentum mz, due to a 
precessing oscillation } about its vertical axis, 
in moving along its path. The lateral com- 
ponent of the ground reaction, which is not 
shown, is simply F = m(r6 + rad), since 
v = rw, while the vertical component is 
approx. N = mg. On taking moments about 
a transverse axis through the center of 
gravity (that is, tangent to its mean path), 
we have + = — Fr, giving 
(A + mr)6 + mrod + Cod = mero. With 
Ad — Cwh = 0, about the vertical, we have 
the oscillation equations for the hoop. 

In the analysis of more complicated systems 
with several degrees of freedom, with the 
direct application of Euler’s equations, it is 
necessary to include the interactions between 
parts of the system, unless a more generalized 
method of analysis, as with Lagrange’s 
equations, is used. Such latter methods are 
not used in this book. 

RUPEN EXKSERGIAN 
The Franklin Institute Laboratories 


ALIPHATIC FLUORINE Compounps, by Alan 
M. Lovelace, William Postelnek and 
Douglas A. Rausch. 370 pages, 6 X 9 in. 
New York, Reinhold Publishing Corp., 
1958. Price, $12.50. 


This book on the preparation and properties 
of aliphatic compounds will be useful to any 
chemist who is interested in the synthesis of 
fluorine-containing organic compounds. 
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The various methods used to introduce 
fluorine into organic molecules are discussed 
concisely and generous reference to the 
literature is evident from the list of 224 
references at the end of this section. Other 
chapters deal with the classical types of 
organic compounds such as alkanes, alkenes, 
alcohols, ethers, etc. Two chapters are 
devoted to nitrogen compounds and these 
contain sections on azomethines, triazines in 
addition to the more common types of com- 
pounds such as amines, nitriles, nitro- 
compounds, etc. A discussion of organome- 
tallic and -metalloid compounds is included. 

Each chapter contains a brief description 
of the various methods used to prepare the 
derivative. A table of compounds, arranged 
according to empirical formula and indicating 
the preparative method by number, the 
yield, physical properties and reference to the 
original literature, is also given in each 
chapter. Altogether these tables contain 
information about 4500 compounds. The 
authors have carefully reviewed the literature 
through 1955 and have included important 
contributions made during the following two 
years, 

Obviously the authors could not give a 
detailed discussion of the properties of all 
aliphatic fluorine compounds but they have 
done a very creditable job of selecting im- 
portant topics for somewhat expanded 
coverage. They have adequately provided 
the reader with references to cover the classes 
of compounds of lesser interest to most 
chemists. 

PAuL TARRANT 
University of Florida 


PROCEEDINGS OF THE SECOND INTERNATIONAL 
CONGREss OF SuRFACE Activity, VoL. I, 
Gas/Liguip AND INTER- 
FACE. 521 pages, diagrams, 6 X 94 in. 
New York, Academic Press Inc., 1958. 
Price: $15.00; complete set of 4 vols., 
$50.00. 


The Proceedings, Volume I, are subdivided 
into the following parts: 1. Insoluble Films; 
2. Soluble Films; 3. Foams; 4. Evaporation 
Retardation; 5. Monolayers Liquid/Liquid 
Interface; 6. Solubilization and Micelles; 7. 
Emulsions; 8. General Phenomena. 

Under these headings, a good review of the 
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modern research on surface activity at the 
International meeting is given, where most 
of the specialists presented their reviews and 
results. The papers that the reviewer con- 
siders most interesting are: H. E. Ries and 
W. A. Kimball on “Structure of Fatty-acid 
Monolayers and a Mechanism for Collapse,” 
in Part 1; B. V. Derjaguin and A. S. Titijev- 
skaya, “Static and Kinetic Stability of Free 
Films and Froths,” in Part 3; Part 4 de- 
scribing the newly developed conservation of 
water by spreading hexadecanol films on the 
surface water reservoirs; the article of M. B. 
Smith and A. E. Alexander on “Solubilization 
of Soap Micelles,”’ in Part 6; and M. Linton 
and K. L. Sutherland on “Dynamic Surface 
Forces Drop Circulation, and Liquid/Liquid 
Mass Transfer,” in Part 8. 

These papers present new results not 
previously published that appear to have 
permanent intrinsic value to the reviewer. 
The whole presentation at the Congress 
presented a very high standard of scientific 
endeavor and all the articles are important 
in their particular field. The whole volume 
is probably the best representation of the 
research in the field of surface phenomena at 
the moment and is well recommended for 
this purpose. 

WLADIMIR PHILIPPOFF 
The Franklin Institute Laboratories 


PROCEEDINGS OF THE SECOND INTERNATIONAL 
ConGrEss oF SuRFACE Activity, II, 
Sottp/Gas INTERFACE, edited by J. H. 
Schulman. 348 pages, illustrations,6 K 94 
in. New York, Academic Press, Inc., 1958. 
Price: $12.60; complete set of 4 vols., 
$50.00. 


This volume on solid-gas interface presents 
the papers of the Congress in two broad 
classifications, physical adsorption containing 
23 papers and chemisorption with 12 papers. 
Although many nations were represented, 
fortunately for the English readers, the con- 
tributions with the exception of a few dis- 
cussions are presented in this language. The 
utility of the collection would have been 
enhanced by having the individual discussions 
immediately following the pertinent papers. 
An index also would be helpful. 

The introductory lecture of Professor G. M. 
Schwab to the section on chemisorption was 
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a most fitting and welcome contribution. 
The proceedings of a Congress of this import- 
ant nature should have such introductions 
giving a brief review of the state of the 
research and a critical summary of the work 
presented, to point out to the young as well 
as other workers in the fields the areas where 
effort might be intensified for greater progress. 
The contributions are given in a very readable 
style. In general, each paper and discussion 
contains an adequate up-to-date set of 
references. While this volume will serve as 
a reference collection of recent contributions 
to the field of research, its value for the 
advanced worker will increase as additions 
are made from future proceedings especially 
if indexes are included. 

R. S. DALTER 
The Franklin Institute Laboratories 


PROCEEDINGS OF THE SECOND INTERNATIONAL 
ConGrEss oF SuRFACE Activity, VoL. III, 
ELECTRICAL PHENOMENA AND SOLID-LIQUID 
INTERFACE, edited by J. H. Schulman. 
621 pages, diagrams, 6 X 94 in. New 
York, Academic Press, Inc., 1958. Price: 
$16.80; complete set of 4 vols., $50.00. 


Volume III of the Proceedings of the Second 
International Congress of Surface Activity 
deals with Electrical Phenomena and the 
Solid-Liquid Interface. There are fifteen 
papers devoted to the first topic and fifty-one 
which deal mainly with the second. The 
author has grouped these later papers under 
six sub-headings, namely: Contact Angles, 
spreading and wetting; Flotation; Adhesion; 
Adsorption ; Lubrication; and Nucleation. 

A large variety of phenomena are presented 
under the heading of electrical phenomena. 
These include studies on the “U-Effect” 
which describe the voltages which arise due 
to mechanical vibrations of the electrical 
double layer, corrosion effects and the kinetics 
of the adsorption of organic materials and 
metal surfaces to name a few. Many of 
these papers are necessarily brief and those 
readers not working directly in these fields 
will find it necessary to do a considerable 
amount of background reading in order to 
appreciate the full impact of many of the 
experiments on physical chemistry theories. 

J. J. Bikerman and J. Gustalla in two 
separate papers point out some oft-times 
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overlooked considerations which surround the 
interpretations of surface tension data. 
These have to do with the apparent in- 
applicability (except in special cases) of the 
Young Equation which relates interfacial 
tension, contact angle and surface tensions 
at solid-liquid interfaces. These papers are 
important since the Young equations (or 
some derivative thereof) are used in inter- 
preting a large amount of current experi- 
mental data. 

The adsorption of low and high molecular 
weight material on metallic and nonmetallic 
surfaces is treated both from the experimental 
and theoretical viewpoints. Despite the 
briefness of many of the articles, they serve 
as informative reviews. 

The effects of surface active media on such 
physical properties as creep, strain and 
rupture are treated in a paper by P.Rehbinder 
and V. Lichtman. This work was performed 
on single metal crystals and thus provides a 
method of correlating adsorption results with 
crystal structure. 

One of the more interesting articles on 
nucleation deals with the formation of ice 
crystals from super-cooled vapors. Thisarticle 
presents recent data obtained on the physics 
of cloud formation and the role played by 
condensation nuclei. 

This volume covers an enormous amount of 
material and the references and the discussion 
sections are most welcome since the phenom- 
ena in general cover both solid state 
physics and physical chemistry. In general 
the articles are well illustrated. This book 
brings together under one cover a group of 
studies of an important class of surface 
sensitive reactions. 

James B. Drew 
The Franklin Institute Laboratories 


PROCEEDINGS OF THE SECOND INTERNATIONAL 
CONGRESS OF SURFACE Activity, VoL. IV, 
INTERFACE (WASHINGS, ETC.) 
AND CELL/WATER INTERFACE, edited by 
J. H. Schulman. 352 pages, diagrams, 
6 X 94 in. New York, Academic Press, 
Inc., 1958. Price: $12.60; complete set of 
4 vols., $50.00. 

As the title implies, the volumes are a 


collection of technical papers presented at 
the Second International Congress of Surface 
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Activity. Vol. IV is subdivided into three 
main sections: Washing, with 17 papers 
followed by a discussion; Suspensions and 
Agglomerations, with 9 papers and a dis- 
cussion; and Cell Water/Interface (Biologi- 
cal), with 13 papers and a discussion. 
Because the volume consists of technical 
papers, the author assumes that the reader is 
working in the same field and there is es- 
sentially no introduction to the subject 
For example, the work on surface activity of 
biological cells appears to be addressed to 
biologists interested in surface activity rather 
than to the whole field of workers who are 
interested in surface activity The editor, 
as he indicated in the Editor’s Note, has 
collected the original papers and presented 
them to the publisher Asa result the papers 
are in three languages; English, French and 
German. Each paper is abstracted in the 
language in which it was written. It would 
have been helpful to present the abstract of 
each paper in the three languages. The 
quality of the papers is excellent with ample 
references to current literature and journal! 
articles. The volume will serve as a reference 
to recent work in the respective fields. 
WILLIAM SCHUELE 
The Franklin Institute Laboratories 


ELEMENTS OF WATER SUPPLY AND WASTE- 
WatTER DtsposaL, by Gordon Maskew 
Fair and John Charles Geyer. 615 pages, 
diagrams, 6 X9 in. New York, John 
Wiley & Sons, Inc., 1958. Price, $8.95. 


Since this book has been written with the 
undergraduate student of sanitary engineering 
in mind, the 175 problems at the end of the 
book illustrating the scientific principles 
presented are an outstanding feature. 

While not as complete as the authors’ 
larger volume on “Water Supply and Waste- 
Water Disposal,’ the present work is entirely 
adequate in presenting the fundamental laws 
involved in the treatment of raw water 
supplies and waste effluents. The two sub- 
jects are often discussed together and much 
duplication is thus avoided. 

Typical of the contents are those chapters 
dealing with “Collection of Surface Water,” 
“Transmission of Water,’ “Distribution of 
Water,” ‘“‘Screening, Sedimentation and Flo- 
tation,” “Filtration,” ‘Biological Treat- 
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ment,” and “Sludge Treatment and Dis- 
posal.’’ As in any textbook of this nature, 
the appendix contains many tables and lists 
of water treatment chemicals that can be 
extremely helpful. A feature of this section 
is a diagram which will be helpful in dealing 
with the Hazen-Williams formula for deter- 
mining flow of water in pipes. 

A. B. MrIpDLETON 
Philadelphia Quartz Co. 


OBSERVATION AND INTERPRETATION, A SyM- 
POSIUM OF PHILOSOPHERS AND PHYSICISTs. 
A Report of the Proceedings of the Ninth 
Symposium of the Colstor Research Society, 
edited by S. Korner. 218 pages, plates, 
74 X 92 in. New York, Academic Press, 
Inc., 1958. Price, $8.00. 


The tape recorder is a boon to the book 
companies. People active in research are 
usually loathe to undertake the job of writing 
a book. However, there is usually no great 
difficulty in the transformation of the tape 
recorder of a conference to print. Unfortun- 
ately, it remains true that the cost in effort 
of a project is, in most cases, closely propor- 
tional to its final value. 

The title, Observation and Interpretation, is 
only weakly descriptive of the symposium 
from which the content of the book under 
review has been taken. The conference 
concerned itself with various philosophical 
problems associated with the interpretation 
of quantum mechanics. The discussions 
dealt mainly with the foundations of the 
concept of probability, the problems of 
objectivity and subjectivity in quantum 
mechanics and the quantum theory of 
measurement. The conference was attended 
by some forty philosophers and physicists. 
Seventeen papers, ranging in length from 
about twenty minutes to about an hour were 
read; these papers constitute somewhat less 
than half of the conference report, the 
remainder being devoted to discussion. 

Although several of the individual papers 
are interesting, the report of the conference 
does not, to my mind, constitute a significant 
book. There is no unifying thread binding 
the discussions of the conference together, no 
general conception provided which would 
allow us to understand why the often con- 
flicting discussions of the conference should 
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find themselves between the same covers. 
The papers were, for the most part, not 
particularly new. The function of the con- 
ference presumably was to permit discussion 
between physicists and philosophers on 
matters of mutual interest. In this, the 
conference was not as successful as one might 
have hoped it would be. 

The book may be of interest to philosophers 
of science in that it reveals something of the 
range of reactions of physicists and philos- 
ophers to quantum mechanics. The book 
also shows that experts can often engage in 
arguments as pointlessly and as fruitlessly as 
amateurs. 

LEONARD EISENBUD 
Bartol Research Foundation 


PROCEEDINGS OF THE COLLOQUIUM SPECTRO- 
SCOPICUM INTERNATIONALE VI, edited by 
W. van Tongeren and F. Freese. 663 

74 X92 in. New 


pages, illustrations, 
York and London, Pergamon Press, Inc., 
1957. Price, $25.00. 


On the 14th through the 19th of May, 1956, 
the Sixth International Colloquium on 
Spectroscopy was held in Amsterdam, Nether- 
lands. This book contains the papers pre- 
sented to that meeting. An adequate review 
of a book of this type is nearly impossible to 
give as the papers are very heterogeneous. 
Experimental work in spectroscopy of the 
ultra-violet, the visible, and the infrared 
regions of the spectra and covering work on 
sources, instruments, and applications are 
presented. Perhaps the fields covered can 
best be shown by listing the titles of the 
sessions as given in the Table of Contents 
and the number of papers under each title: 
(1) Plain Spectroscopy (10 papers); (2) 
Emission Spectroscopy (8 papers); (3) Arc 
Emission (17 papers); (4) Spark Emission 
(15 papers); (5) Emission (General) (9 
papers); (6) Instrumentation (Emission) (19 
papers); (7) Instrumentation of Absorption 
Methods (8 papers) ; (8) Ultra-Violet Absorp- 
tiometry (7 papers); and (9) Infrared 
Spectroscopy (15 papers). 

These papers are presented in French, 
German and English by scientists from twenty 
countries including Western Europe, Russia 
and the satellites, China, South Africa, and 
the U. S. A. 
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Following each paper, questions and com- 
ments from the floor are recorded. Many of 
the papers are abstracted in English and some 
in English, French and German. There are 
about a dozen papers in which only the 
summaries are presented, as the papers were 
not received by the editors. 

This volume contains such a variety of 
papers presenting new techniques for all the 
fields of spectroscopy that a study of this 
work should be a necessity to all spectrog- 
raphers and designers of spectrographs. 

LuTHER C. Conant, JR. 
The Franklin Institute Laboratories 


Notes ON ANALOG-DIGITAL CONVERSION 
TECHNIQUES, edited by J. W. Brean, M. S. 
Osman, D. T. Ross, J. E. Ward, and A. K. 
Susskind. 410 pages, diagrams, 53 X 9 in. 
New York, John Wiley & Sons, Inc., 1958. 
Price, $10.00. 


Five members of M.1.T.’s Servomechanisms 
Laboratory have combined efforts to treat 
this specialized subject. Perhaps ‘“‘special- 
ized” is not the correct word because the 
techniques of analog-to-digital conversion now 
cut across many technical fields. These in- 
clude such widely separated practices as 
radar, computation, and numerical control of 
machinery. The analog quantities treated 
may be voltage signals as well as physical 
motions, and in the two cases the means of 
conversion to digital form are divergent 
indeed. 

Casual use of the term “digital” is also an 
oversimplification. What is implied is the 
sampling and quantizing of the information 
with the eventual application of some form 
of code. It can be seen that the theoretical 
background of sampling and the principles of 
code development are absorbing subjects in 
themselves. These are taken up by the 
authors before turning to the detailed con- 
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siderations of analog-to-digital conversion 
mechanisms. 

The engineering considerations confronting 
the designer in the instance of a numerically- 
controlled machine tool are intricate but 
interesting. Questions of electronic circuits, 
optics, magnetics, and mechanics of extra- 
ordinary accuracy may be involved. The 
book attempts to give a_ state-of-the-art 
picture of each specific technique involved. 
Decoding of the digitally stored information 
or instructions is discussed following the 
means of coding. Finally the over-all design 
of digital instrumentation systems is taken 
up to show how the components may affect 
each other. 

The book appears as an informal publica- 
tion employing offset printing and showing 
some lack of good editing. Parts are wordy, 
as though ideas and information were hard 
to find, while some portions are brief and 
unclear. These comments do not, however, 
reflect upon the authors’ technical com- 
petence, of which there is ample evidence 
through the points raised in the text. Rather, 
there is the suggestion that the whole publica- 
tion project was hurried and had to depend 
upon the spare time busy men could snatch 
from their other work. One serious short- 
coming, particularly in a book like this, 
where terms and definitions are so important, 
is the absence of any kind of index. However, 
the reviewer has already found the text and 
illustrations useful in his work, and it is 
certainly to be expected that others will. It 
might be pointed out also in closing that this 
is the group that did the first serious work on 
the numerical control of large milling ma- 
chines. The very extensive projects carried 
out in the Servomechanisms Laboratory 
undoubtedly provide much of the experience 
that is assembled in this impressively compre- 
hensive book. 

C. W. HARGENS 
The Franklin Institute Laboratories 
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SCIENTIFIC PROGRAMMING IN BUSINESS AND 
InpustRy, by Andrew Vazsonyi. 474 
pages, diagrams, 6 X 9 in. New York, 
John Wiley & Sons, Inc., 1958. Price, 
$13.50. 


Written for the businessmen who want to 
use scientific techniques in solving managerial 
problems, this book sets forth in understand- 
able language, the ways in which mathematics 
can aid management—linear programming, 
statistical decision theory and operations 
research. The author carefully presents 
fundamentals and then goes into detail about 


their use in practical business problems such . 


as transportation allocation, production, 


inventory, warehousing, stable employment, 
assembly line scheduling, accounting, etc. 
In these days, no business executive can afford 
to ignore the techniques of mathematics now 
being applied so successfully in all lines of 


business ; if he knows nothing about them, he 
should get this book and study it carefully. 


VEcToR SPACES AND MATRICES, by Robert M. 
Thrall and Leonard Tornheim. 318 pages, 
6 X 9in. New York, John Wiley & Sons, 
Inc., 1957. Price, $6.75. 


In an attempt to bridge the gap between 
beginning and higher levels of the theory of 
vector spaces and matrices, the authors have 
used a dual approach on two levels—one 
concrete (via matrices) and the other axio- 
matic (via linear transformation). Eleven 
chapters cover the field thoroughly, including 
vector spaces, linear transformations, deter- 
minants, orthogonal and unitary equivalence, 
polynomial rings, linear inequalities, etc. 
The theory of games is briefly touched on in 
the final chapter. 


Circuit ANALYSIS OF TRANSMISSION LINEs, 
by John Stewart. 186 pages, diagrams, 
6 X 9in. New York, John Wiley & Sons, 
Inc., 1958. Price, $5.50. 


This new addition to the literature provides 
a compact but comprehensive treatment of 
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the ordinary transmission line, by the usual 
approximation to a “‘lossless’’ line. The 
Maxwell equations are not gone into, hence 
the test is suitable for instruction at the 
undergraduate level. It is well arranged for 
this purpose, taking up in turn the derivation 
of the basic equations for “‘lossless’’ and 
“lossy”’ lines, wave propagation, resonant 
lines, and cavities, equivalent lumped circuits, 
standing-wave ratio and impedance measure 
ments, and the use of the Smith chart for 
graphical solutions. Basic data for trans- 
mission line calculations are included in a 
five page appendix. While the emphasis is 
on radio-frequency lines, the audio and power 
engineer will also find the book of interest 
and assistance. 


H1GH-SpEED Data ProcessinG, by C. C. Got- 
lieb and J. N. P. Hume. 338 pages, dia- 
grams, 6 X 9in. New York, McGraw-Hill 
Book Co., Inc., 1958. Price, $9.50. 


The authors have thoroughly outlined the 
general techniques of processing data at high 
speeds, particularly for business purposes. 
This text covers a wide range of subjects, 
from the method of representing information 
in a processor to advances in automatic 
programming. Not limiting itself to any 
particular machine, the book covers data 
processing of all types and from all manu- 
facturers, using a hypothetical machine as a 
working model. 


MoRPHOGENESIS: THE ANALYSIS OF MOL- 
LUSCAN DEVELOPMENT, by Chr. P. Raven. 
311 pages, diagrams, 53 X 8} in. New 
York, London, Paris, Pergamon Press, Inc., 
1958. Price, $10.00. 


This volume is a survey, both descriptive 
and experimental, of the present state of 
knowledge on development in the molluscs. 
The chapters deal mainly with the causal 
analysis of morphogenesis. Cytochemical 
and experimental embryological investiga- 
tions on molluscs have been completely 
covered, and cytological aspects of develop- 
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ment are especially emphasized. Researches 
from the author’s own laboratory claim an 
important place in this book. 


HicH Qua.ity Sounp REPRODUCTION, by 
James Moir. 591 pages, plates, diagrams, 
54 X 84 in. New York, The Macmillan 
Company, 1958. Price, $14.00. 


In this profusely illustrated volume on 
sound reproduction, the author discusses the 
different characteristics of acoustic conditions 
of auditoria. Considerable space is devoted 
to the design of amplifiers and to a discussion 
of loudspeaker mountings and enclosures. 
An appendix is included at the end of each 
chapter for easy reference. The professional 
engineer, as well as the amateur interested in 
the problem of reproducing sound with a high 
degree of fidelity, will find this book provides 
a wealth of material. 


Puysics FOR ENGINEERS AND SCIENTISTS, by 
Richard G. Fowler and Donald I. Meyer. 
546 pages, illustrations, 6 X 9in. Boston, 
Allyn and Bacon, Inc., 1958. Price, $8.00. 


A “new look” physics text which discards 
the classical division of subject matter for a 
unified approach and emphasizes modern 
insights and developments. Calculus is used, 
but may be omitted for more elementary 
courses. Instructors should appreciate the 
convenient division into 30 topical chapters 
averaging 17 pages and suitable for a week's 
work each. Students will appreciate the 
wealth of excellent illustrations. 


TABLE FOR THE SOLUTION OF CuBIC Equa- 
TIoNs, by Herbert E. Salzer, Charles H. 
Richards and Isabelle Arsham. 161 pages, 
54 X 8in. New York, McGraw-Hill Book 
Co., Inc., 1958. Price, $7.50. 


This useful table makes it possible to find 
the three roots for any cubic equation having 
real coefficients, in a few minutes with a desk 
calculator. More accurate and easier to use 
than other existing tables, this table will be a 
valuable aid for engineers, physicists and 
applied mathematicians. Adequate instruc- 
tions for using the table are given in the intro- 
duction, including interpolation. The table 
covers every possible set of real coefficients. 
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AGRICOLA ON METALS, edited by Bern Dibner. 
128 pages, illustrations, 84 X 103 in. 
Conn., Burndy Library, 1958. Price, $2.50. 


The mining and metallurgical industries of 
16th century Germany is the subject of 
Georgius Agricola's book. This historic 
work, originally titled “De Re Metallica” 
and written in Latin, was translated by 
Herbert Hoover and his wife in 1912 after 
five years of research. Reproducing over 70 
selected illustrations from the original edition, 
the Editor of the present offering, Bern 
Dibner, presents a distillation for the reader 
interested primarily in the historical signifi- 
cance of the origin of metals. This text is 
a record of a highly developed industrial 
culture that flourished some 200 years ahead 
of what we call ‘the industrial revolution.” 


PRINCIPLES OF GEOCHEMISTRY, by Brian 
Mason. Second edition, 310 pages, dia- 
grams, 6 X 9 in. New York, John Wiley 
& Sons, Inc., 1958. Price, $8.50. 


A pleasantly written systematic account of 
the chief facts and ideas of the chemistry of 
the earth and of its chemical and physical 
evolution, adapted from notes for a one- 
semester course taught by the author, who is 
Curator of Physical Geology and Mineralogy 
at the American Museum of Natural History 
and a Professor at Columbia. It should be 
useful to students of geology, astronomy, and 
chemistry. There is a substantial bibliogra- 
phy, with annotations. 


FUNDAMENTALS OF CHEMICAL ENGINEERING, 
by Maurice G. Larian. 644 pages, dia- 
grams, 6 X 9in. Englewood Cliffs (N. J.). 
Prentice-Hall, Inc., 1958. Price, $10.50. 


A new text for undergraduates, this book 
emphasizes the analytical rather than the 
empirical approach. Three main sections 
cover General Principles (flow of fluids, heat 
and mass transfer, evaporation), Separation 
of Mixtures by Interphase Mass Transfer 
Operations (drying of solids, distillation 
techniques, extraction and gas absorption), 
and Mechanical Separation of Heterogeneous 
Mixtures (gravitational and centrifugal 
methods and filtration). The author, Pro- 
fessor of Chemical Engineering at Michigan 
State University, has been guided by his 
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own classroom experiences in the preparation 
of the text. 


Atomic ENERGY LEVELS, As DERIVED FROM 
THE ANALyYsis OF OpTicaL SpecTRA, VOL. 
III, by Charlotte E. Moore. 245 pages, 
9 X 114 in. Washington, U. S. Govern- 
ment Printing Office, 1958. Price, $2.50; 
to foreign countries, $3.13 (in U. S. 
exchange). 


Volume III in this series covers the atomic 
energy levels of the elements Molybdenum 
through Lanthanum (Z = 42 through 57) 
and Hafnium through Actinium (Z = 72 
through 89). It includes 124 spectra, with 
the data presented as in Vols. I and II. 
With the publication of the present volume, 
the tabulation of the Periodic Table is com- 
plete except for two groups of rare-earth 
elements which will be published as Vol. IV. 
Until this last volume is available, a current 
list of references to analyses of rare-earth 
spectra is included in Vol. III. The series 
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is published by the National Bureau of 
Standards. 


SERIAL PUBLICATIONS OF THE SOVIET UNION 
1939-1957, compiled by Rudolf Smits, 
459 pages, 9 X 114in. Washington, U. S. 
Government Printing Office, 1958. Price, 
$2.75. 


This bibliographic checklist of Soviet serial 
publications includes government and non- 
official publications appearing either regularly 
or irregularly in the Soviet Union in 1939 or 
later—except those in oriental languages. 
The checklist, published by The Library of 
Congress, includes the title, place of publica- 
tion, issuing body and year the publication 
began, as well as the names of U. S. libraries 
holding the publication. A subject guide 
referring to specific publications will help the 
English reader immeasurably. With the 
increasing emphasis on Soviet work, this 
checklist will be a valuable aid to scientists 
in this country. 
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MATHEMATICAL TABLES AND FORMULAE, by F.J.Camm. 144 pages, 44 X 7}in. New York, 
Philosophical Library, 1958. Price, $2.75. 

OPERATIONAL MATHEMATICS, by Ruel V. Churchill. Second edition, 337 pages, diagrams, 
6 X9in. New York, McGraw-Hill Book Co., Inc., 1958. Price, $7.00. 

DIFFERENTIAL Equations, by Forest Ray Moulton. Unabridged and unaltered republication 
of the first edition, 395 pages, 5 X 8 in. New York, Dover Publications, 1958. Price, 
$2.00 (paper). 

MATHEMATICAL TABLES, by Herbert Bristol Dwight. Second edition, 220 pages, 5} X 8 in. 
New York, Dover Publications, 1958. Price, $1.75 (paper). 

MICROWAVE TRANSMISSION DEsIGN Data, by Theodore Moreno. Unabridged and unaltered 
republication of the first edition, 248 pages, 5 X 8 in. New York, Dover Publications, 
1958. Price, $1.50 (paper). 

History oF Matuematics, by D. E. Smith. Unaltered and unabridged republication of the 
last edition. New York, Dover Publications, 1958. Vol. I, 596 pages, 54 X 8 in. Vol. 
II, 725 pages. Price, $2.75 each or $5.00 for both. 

ConcIsE INTERNATIONAL DICTIONARY OF MECHANICS AND GEOLOGY (English-French-German- 
Spanish), by S. A. Cooper. 400 pages, 43 X 7} in. New York, Philosophical Library, 
1958. Price, $6.00. 

Aupio DesiGN HANDBOOK, by H. A. Hartley. 224 pages, diagrams, 54 X 84in. New York, 
Gernsback Library, Inc., 1958. Price: $2.90 (paper) ; $5.00 (hard cover). 

ELECTRONIC MEASURING INSTRUMENTS, by E. H. W. Banner. Second edition, revised, 
496 pages, diagrams, 54 X 8} in. New York, The Macmillan Company, 1958. Price, 
$7.95. 

INTRODUCTION TO NUCLEAR ENGINEERING, by Richard Stephenson. Second edition, 491 pages, 
illustrations, 6 X 9 in. New York, McGraw-Hill Book Co., Inc., 1958. Price, $9.50. 

CuHemistry, by James V. Quagliano. 877 pages, diagrams, 6 X 9 in. Englewood Cliffs 
(N. J.), Prentice-Hall, Inc., 1958. Price, $6.95. 

A COMPREHENSIVE BIBLIOGRAPHY ON OPERATIONS RESEARCH. 188 pages, 8$ X 11 in. 
New York, John Wiley & Sons, Inc., 1958. Price, $6.50. 

A SELECTED BIBLIOGRAPHY OF RESEARCH AND DEVELOPMENT AND ITS IMPACT ON THE ECONOMY. 
21 pages, 7? X 10} in. Washington, D. C., National Science Foundation, 1958. Price, 
not given. 

TABLES OF TEMPERATURE, RELATIVE HUMIDITY AND PRECIPITATION FOR THE WORLD. 
6 X 92 in. London, Her Majesty’s Stationery Office, 1958. One volume (composite), 
or six parts, covering North America, South America, Europe, Africa, Asia and Australasia. 


Various paging, various prices, paper. 
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Space Technology Tools.—New 
space technology tools, developed to 
help produce vehicles for flight into 
outer space, have been disclosed by 
the General Electric Company. In- 
cluded in the family of new G-E 
space technology tools is America’s 
biggest plasma jet generator for 
testing space models under the high- 
temperature and the chemical condi- 
tions of atmospheric re-entry from 
outer space, a unique, 120-ft. long 
shock tunnel, largest in the country, 
used to study space vehicle design 
requirements at velocities up to 25 
times greater than the speed of sound, 
a solar furnace that automatically 
tracks the sun to provide clean, 
radiant heat energy for making non- 
contaminating tests, and electronic 
mass-accelerators—which will hurl a 
mass at velocities over twice as great 
as those produced by a 16-in. cannon 
—for free-flight tests. 

According to George F. Metcalf, 
top General Electric space flight 
executive and general manager of 
the Company’s Missile and Ordnance 
Systems Department, ‘these elec- 
tronic, atomic, high-temperature, and 
high-speed probing devices support 
the modern space pioneer in the same 
way the axe and flintlock served the 
early American frontiersman.”’ 

Located at the Company’s Missile 
Ordnance Systems Department, in 
Philadelphia, some of the equipment 
is being used to help solve the nose 
cone re-entry problem on the Air 
Force’s ATLAS and THOR ballistic 
missile programs. The equipment is 
operated by engineers and scientists 
in G-E’s Aerosciences Laboratory. 
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Newest of the space technology 
tools is the giant plasma generator, 
just now completing its development 
stage. Installed at G-E’s Switchgear 
Department, one of the few facilities 
in the free world capable of generating 
the electric current it needs, the 
plasma jet yields temperatures more 
than twice the surface temperature 
of the sun. 

The large, air-stabilized electric arc, 
enclosed by a chamber 18 in. through 
and 3 ft. high, as contrasted to fore- 
runners of 4 by 10 in., is a three-phase 
a-c. model and has been run at the 
rate of 15,000 kw. At this rate of 
usage, based upon the average daily 
home use, the power needed to operate 
the plasma generator at its peak is 
equivalent to that flowing into 5,400, 
000 homes. 

“‘Advantages of an arc of this size 
are that more nearly full-scale models 
of re-entry vehicles can be tested in 
its plasma flow and that the simulated 
Mach 12-25 speed of the plasma 
jet more nearly approximates actual 
space-flight conditions,’’ Metcalf said. 

Much smaller stabilized arcs, an- 
nounced by G-E over a year ago, 
have been used for testing materials 
in small-scale models and have as- 
sisted engineers to study ionic propul- 
sion methods and upper atmospheric 
conditions. G-E, a pioneer in stabi- 
lized arc development, has designed 
and built both liquid and gas stabi- 
lized arcs. 

Stabilized arcs, whether liquid or 
gas, differ from normal electric arcs 
in that the arc column is controlled, 
or stabilized, by a squeezing blanket 
of liquid or gas. This blanket, forced 
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into the arc container tangentially, 
develops a hollow center, or vortex, 
through which the arc is struck. The 
arc column is restricted by the size of 
the vortex, and increased current can 
be fed into the arc, raising the tem- 
perature correspondingly. A plasma 
jet is then forced out to test materials 
at high speeds and high temperatures. 
In conventional arcs, increased cur- 
rent causes the arc column to broaden, 
spreading the heat over a wider area; 
hence the value of stabilization for 
high heat concentration. 

Perhaps the most impressive of 
G-E’s space technology tools is this 
country's largest shock tunnel, 6 in. 
by 120 ft., which permits high- 
temperature and high speed studies— 
aerothermodynamics. Unlike other, 
conventional shock tubes, this space 
technology tunnel is equipped with 
an expansion nozzle and a large 
reservoir, which contains the speci- 
men. Studied by means of Schlieren 
photography, specimens are subjected 
to Mach 15-25 blasts of air that are 
compressed and driven, by exploding 
gases, down the tube. The 18,000° F. 
and 5000 psi. air then flows through 
the expansion nozzle in front of the 
reservoir over the specimen. Flow 
patterns set up by these fast-moving 
gases yield data that will largely 
determine spacecraft design. 

Three smaller, companion tubes, of 
varying lengths, diameters, and inside 
configurations, support G-E’s large, 
space-age shock tunnel by eliminating 
smaller testing problems and by 
providing basic information on high- 
temperature, gaseous physics. 

For studies requiring relatively low 
temperatures, aeroscientists at G-E 
use a solar furnace. This radiant 
energy source has the advantage of 
providing clean heat cheaply. It 
takes heat directly from the sun, or, 
on cloudy days, from a large arc light. 
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When using the sun as a source of 
heat, the heliostat-equipped solar 
furnace is fed by the sun’s rays 
automatically ; for the heliostat, using 
mirrors, directs the sun into the 
furnace’s concave reflector, where its 
3000° C. heat is focused onto a speci- 
men. A hollow quartz hemisphere, 
capable of being evacuated or of 
allowing gas mixtures that simulate 
space vehicle environments, further 
enhances the value of this otherwise 
simple mechanism. 

Other members of the family of 
space technology tools have been con- 
ceived. An example is the arc dis- 
charge hypersonic gun, operated by 
striking arcs between evenly spaced 
electrodes along a tube to heat a 
gas—helium, for example—whose ex- 
pansion will force a projectile at 
increasing speeds down the tube. 
The uniformly increasing pressure 
behind the mass should propel it, by 
the time it reaches the end of the 
tube, at a muzzle velocity vastly 
greater than that produced by a 
16-in. cannon. 

Chemical means of propulsion, for 
example, yield less than 25,000 fps., 
while electrical means theoretically 
may push a mass to 60,000 fps. The 
purpose for such high velocities is to 
provide free flight tests after the mass 
leaves the tube at speeds equivalent 
to those in space. The mass so pro- 
pelled, possibly a model space vehicle, 
can subsequently be studied, by 
spectroscopic or telemetric means, to 
provide information on actual space 
conditions. 


Device Measures’ Transistor 
Switching Time.—A device for testing 
transistor switching times incorpor- 
ating a new precision technique for 
measuring durations in the millimicro- 
second region has been announced by 
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Atronic Products, Inc., Bala-Cynwyd, 
Pa. 

Designed specifically for analyzing 
the time characteristics of high speed 
transistors, the API 300 will measure 
the switching time of the leading or 
trailing edge of the output pulse with 
an error of 3.0 millimicroseconds in 
the 3 to 100 millimicrosecond range. 
Direct meter reading featuring a pre- 
settable go-no-go level makes the 
unit highly adaptable to production 
testing as well as to laboratory 
application. 

In operation, the transistor under 
test is inserted into the input connec- 
tor of the API 300. Self contained 
transistor bias voltages are provided, 
and an extremely precise driving 
pulse is applied to the transistor. 
By sampling that portion of the out- 
put pulse to be measured with a pulse 
whose characteristics are known, time 
is essentially converted to a voltage 
whose amplitude is proportional to 
duration. This voltage is amplified 
and displayed on a meter in terms of 
millimicroseconds. Functionally, the 
API 300 consists of : 

1. Input Section which holds the 
transistor and makes the proper input 
and output connections. 

2. Power Supply which provides all 
voltages including the transistor bias 
voltages. The unit operates from a 
60 cycle, 110 volt source. 

3. Pulse Generator. This section 
generates a precision pulse with a 
rise time of 1.0 millimicrosecond. 
This pulse is used to drive the transis- 
tor. It is also fed into the Delay 
Section. 

4. Delay Section consists of an 
extremely accurate delay line which 
delays the precision pulse a fixed 
amount before being fed into the 
Gate Section. 

5. Gate Section is triggered by the 
pulse from the Delay Section. It 
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samples the transistor output; inte- 
grates it and produces a d-c. level 
proportional to the switching time. 
This signal is amplified and fed to the 
output section. 

6. Output Section. The output of 
the Gate Section is a voltage level 
whose amplitude is proportional to 
time. This is fed into a voltmeter 
whose dial is calibrated in millimicro- 
seconds for the direct reading of the 
durations to be measured. By this 
technique, extremely precise measure- 
ments can be made in the 3 to 100 
millimicrosecond region. 

The model API 300 is available 
either as a bench top unit or in a rack 
mounting version. Other models are 
available which will separate carrier 
storage time from fall time and pro- 
vide measurements in different ranges 
of switching speeds. 


New Dust Control Equipment and 
Chem-Jet Systems.—New equipment 
especially developed for the control 
of fine dust will be exhibited by The 
Johnson-March Corp., Philadelphia 
dust control engineers, at the Ameri- 
can Mining Congress Show to be held 
in San Francisco from September 22 
through 25. 

Featured at the dust control exhibit 
will be new equipment for the control 
of dust at discharge of dry collectors 
and for fly ash conditioning, sintering 
discharge, bin unloading, rail car 
unloading, and similar applications. 

A working model of a Johnson- 
March unit known asa Type A Hydro- 
Precipitator Scrubber will be shown. 
This equipment is used for the scrub- 
bing of gases containing dusts ranging 
below five microns in particle size and 
for fume control. It is specifically 
designed for use with roasters, dryers, 
pulverizers, and in smelting and 
sintering plant operations. 

The scrubber is said to perform at 
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99.8 per cent efficiency with dust 
loadings at five grains per cubic foot 
and with 70 per cent of the material 
below the five-micron-size. 

The company’s Chem-Jet System 
for controlling dust in metallic and 
non-metallic mining operations also 
will be on display. This unit sup- 
presses dust at its source before it 
can become airborne instead of col- 
lecting the dust as older mechanical 
or dry systems do. Multiple jets 
apply very small quantities of a sur- 
face-active solution containing Com- 
pound M-R, a compound specially 
formulated to reduce surface tension. 
An average of 1} gallons of solution 
per ton of material is added. This 
4 of 1 per cent added moisture does 
not actually wet or adversely affect 
the raw material treated. 

Also on exhibit will be Compound 
SP, a solution for outdoor storage 
protection. Samples taken from piles 
of various raw materials which were 
treated by the compound will be 
shown. These samples will show the 
binding effect of the compound on the 
material, shallow depth of penetration 
of the compound, and unweathered 
condition of the protected material. 
Compound SP does not affect later 
use of the stored material. 


Reactor Safety Fuse.—A reactor 
safety ‘‘fuse,’”’ which automatically 
shuts down high power nuclear re- 
search reactors if operating control 
is lost, has been designed, built, and 
‘successfully tested for the Atomic 
Energy Commission by Atomics In- 
ternational, a division of North Ameri- 
can Aviation, Inc. 

The fuse is a safety device designed 
to provide absolute protection by 
shutting down ‘‘pool’’ type research 
reactors without the use of external 
controls should an abnormal operating 
condition occur. The development is 
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part of the Commission’s reactor 
safety program. 

The fuse was tested recently at the 
Commission’s experimental facility, 
known as Special Power Excursion 
Reactor Test (SPERT), at the Na- 
tional Reactor Testing Station in 
Idaho. It was subjected to severe 
power surges or  ‘“‘transients.”’ 
SPERT, operated by Phillips Petrol- 
eum Company, was designed to study 
the characteristics of power surges in 
nuclear reactors. 

The fuse prevented the reactor 
power from rising excessively during 
a severe transient in which the rate 
of power increase was tripling each 
one-hundredth of a second. The fuse 
terminated the power surge at a level 
25 times lower than the reactor 
would have reached without the 
device. 

In the test, the total energy gener- 
ated by the reactor was less than 
one-tenth that which would have 
been produced without the fuse. 
Energy release is an indication of the 
potential hazard resulting from an 
abnormal reactor operating condition. 

The fuse is completely self-con- 
tained and requires no operator or 
external connections. It includes a 
cylinder, part of which is in the 
reactor core. The part of the cylinder 
outside the core contains boron tri- 
fluoride gas. In the event of a 
reactor transient, or an excessively 
high level of operation, the gas is 
automatically released into the other 
part of the cylinder. The gas acts 
as a reactor “poison” by absorbing 
neutrons needed to maintain the nu- 
clear chain reaction. The gas transfer 
brings the reactor under control 
almost instantaneously and without 
damage to the reactor. 


Medical Photography.—A process 
which provides better medical pho- 
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tography in capillary structure se- 
quence pictures—through use of an 
instrument adapted from another 
scientific field—has been developed 
at the University of Wisconsin 
Hospitals. 

Greater clarity and illumination in 
capillary photos were achieved by 
using a metallurgic microscope as 
source of light. The new process 
replaces the old method of a quartz 
rod as light supply in the capillaro- 
scope. The combination was de- 
veloped by Dr. Enrique Valdivia, 
instructor in pathology at the UW 
Medical School, with the technical 
assistance of Homer Montague, de- 
partment of photography director at 
the hospitals. 

The new process provides a full 
diffusion of light without heat while 
the old method of focused light 
through a quartz rod tip limited the 


lighted area and produced intensive 
heat in the tip which required con- 
tinuous liquid cooling. The need for 


cooling is eliminated in the new 
procedure. 

The innovation provides sharp de- 
tail in the pulmonary structure and 
permits greater enlargement without 
loss of clarity. 

The process was developed for use 
in recording the capillary structure 
change in the lung of a guinea pig 
under simulated high altitude condi- 
tions, part of a continuing study of 
high altitude acclimatization at the 
UW Medical School. 


Combination Immunization Shot.— 
A four-in-one shot designed to im- 
munize children against poliomyelitis, 
whooping cough, diphtheria and 
tetanus has been tested successfully 
in approximately 300 Detroit children, 
according to an article in the June 28 
issue of the Journal of the American 
Medical Association, 
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Children who participated in the 
study ranged in age from 2} months 
to 5 years. Seventy-five per cent, or 
224 of them, completed the series of 
three monthly inoculations, and com- 
parisons of bloods obtained before and 
after the inoculations showed favor- 
able response, the article says. 

The product, known by the trade- 
mark Quadrigen, is not yet available 
for commercial distribution. Devel- 
oped by Parke, Davis & Company, it 
is the result of two years’ research on 
combining polio vaccine with a three- 
in-one shot, diphtheria and tetanus 
toxoids and whooping cough vaccine 
(DPT), now in general use by 
pediatricians. 

The purpose of the research pro- 
gram was to develop a vaccine which 
would protect children against polio, 
as well as the other 3 diseases, with- 
out increasing the volume of vaccine 
or number of injections the children 
now receive for diphtheria, tetanus 
and whooping cough. 

By a special series of manufacturing 
processes, the four components are 
concentrated and provide effective- 
ness with three doses of 3 cc. (just 
over 8 drops) of Quadrigen for each 
inoculation. The Detroit trial was 
conducted to test the response of 
children to the combined vaccine. 

A significant finding in the lengthy 
study was the low incidence of 
naturally acquired resistance to polio 
in 7-9 month old children, indicating 
a need for vaccine inoculation prior 
to this age. Injections of Quadrigen 
produced good antibody response to 
the polio component, as well as to 
the other 3 antigens of the four-way 
vaccine in children from 2} months 
to five years of age. A fourth dose 
was found to be desirable for children 
younger than four months. 

Three Detroit Health Department 
physicians—Dr. Joseph G. Molner, 
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Health Commissioner; Dr. C. D. 
Barrett, Jr., Director of Maternal and 
Child Health ; and Dr. C. P. Anderson, 
Deputy Commissioner of Health, 
Medical Service—conducted the trial 
in conjunction with Dr. E. A. Timm, 
Dr. B. I. Wilner, Dr. H. E. Carnes 
and Dr. I. W. McLean, Jr., members 
of the Parke-Davis research staff. 
The seven doctors collaborated in 
writing this first report on Quadrigen 
for the Journal of the American 
Medical Association. 

Children who participated in the 
study are currently receiving a routine 
“‘booster’’ dose of Quadrigen 15 to 16 
months following their primary series 
of inoculations. Antibody evaluation 
studies of this booster will be a part 
of a more complete report at a later 
date, according to the authors. 


High-Vacuum Thermionic Conver- 
ter.—A working model of a thermionic 
converter—a device that produces 
electricity directly from heat—has 
been described by scientists of the 
General Electric Research Laboratory, 
Schenectady, N. Y. The develop- 
ment was revealed at a meeting of the 
American Physical Society at Cornell 
University. 

The new device, a combination of 
metal and ceramic disks surrounding 
a high vacuum, is the size of a quarter, 
and produces electricity when the 
flame of a blowtorch is played upon it. 
Last year the GE Research Labora- 
tory announced the results of Dr. 
Volney C. Wilson's original work on 
the thermionic converter. The new 
converter, with its high vacuum, is in 
contrast to the original, which was 
filled with gas. 

“In the relatively short time that 
has elapsed since the announcement 
of the original converter, we have 
made great progress in exploring the 
capacities and limitations of such 
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devices,’’ declared Dr. Virgil L. Stout, 
manager of the Physical Electronics 
Section of the Research Laboratory. 
‘“‘In the process we have gone from 
an experimental laboratory device to 
a working model, and there is no 
doubt that further research and 
development will lead us not only to 
more complete scientific understand- 
ing, but also to practical applications.”’ 

Possible applications, he added, 
could be found wherever a high tem- 
perature source of heat is available 
and a supply of electricity is needed. 
Any source of heat, nuclear or con- 
ventional, could be used. The con- 
verter, as a result, should prove 
valuable in a wide variety of applica- 
tions, including aircraft, missiles, and 
satellites. It is estimated that con- 
verters of this size will be made 
capable of operating in the one- to 
ten-watt range. 

The new converter is constructed 
by means of techniques similar to 
those used for high temperature 
electronic tubes, according to a paper 
co-authored by Dr. Harold F. Webster 
and James E. Beggs. Beggs, the 
inventor of the device, worked closely 
with Webster, who performed theo- 
retical analyses. This type of design, | 
Webster pointed out in presenting the 
paper, has the advantage of operating 
at a lower temperature than the gas- 
filled converter (1500° F. vs 3000° F.), 
which makes material requirements 
easier to meet. 

In the thermionic converter, two 
electrodes are held at high, but 
different, temperatures. Electrons 
are ‘‘boiled out” of the hotter cathode 
and are collected by the relatively 
cool anode. They can then flow 
through an external circuit and do 
work. 

The original announcement of Wil- 
son’s converter referred to ‘new 
approaches” to the gas environment 
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within the foot-long glass tube. Dr. 
Wilson revealed at the Cornell meet- 
ing that the gas involved is cesium 
vapor, and that the vapor partially 
ionizes, thus neutralizing the space- 
charge effect that would otherwise 
block the flow of electrons from 
cathode to anode. 

Equally important, he added, is 
the fact that the cesium coats the 
electrodes, thus affecting their work 
functions (the amount of energy that 
must be used in order to free an 
electron from the electrode and release 
it into the surrounding space). The 
difference between the work functions 
of the cathode and the anode is the 
key to the successful operation of the 
thermionic converter, Wilson revealed 
since this difference provides the 
voltage that drives the electrons 
through the circuit. 

The problem of reducing the in- 
ternal resistance of the converter 
can be approached in several ways, 
he pointed out. One is by using a 
partly ionized gas, as he did in his 
original model. A second is by 
placing the electrodes extremely close 
together, which is the method used 
in the new high vacuum converter. 


Standard Case Improves Pulse 
Transformer Packages.—A new 
standard plastic case specifically de- 
signed to eliminate the mounting and 
encapsulating problems of a wide 
range of Pulse Transformers has been 
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announced by the Milton Ross Metals 
Company, of Southampton, Pa. The 
case, ajcomplete pulse transformer 
package rather than just the encap- 
sulation] shell of former types, has 
four special features: Custom fit for 
all sizes of transformers ; unobstructed 
access for assembling and wiring the 
components; a leak resistant press-fit 


seal; and a _ special keyway for 
identifying connectors and inserting 
the package. 


Available in a variety of military 
approved materials which permit oper- 
ations at temperatures up to 500° F., 
the new case is composed of two 
matched parts: A circular shell for 
holding the potting compound and a 
flat plastic header with connecting 
pins molded-in for mounting and 
wiring the component. 

The header, or PMS (Plastic to 
Metal Seal), is flanged to fit a recessed 
ring in the shell. This provides a 
three-surface, press-fit seal which min- 
imizes leakage of the encapsulating 
material. For certain types of trans- 
formers the PMS contains a threaded 
metal insert for the cup-core mounting 
screw. 

To provide custom fit for a wide 
range of transformers the cases are 
available in 25 different diameters 
ranging from 3; in. to 24 in. There 
is also a choice of 25 different heights 
from 3% in. to 24 in. Cases may be 
specified with 4-5-6-7-8-9 or 10 pins 
in various pitch-circle diameters up 
to 13 in. 
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